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ABSTRACT 


A postflight FORTRAN program called “radar” reads and analyzes ground-based radar data. The output includes 
position, velocity, and acceleration parameters. Airdata parameters are also provided if atmospheric characteristics 
are input. This program can read data from any radar in three formats. Geocentric Cartesian position can also be 
used as input, which may be from an inertial navigation or Global Positioning System. Options include spike 
removal, data filtering, and atmospheric refraction corrections. Atmospheric refraction can be corrected using the 
quick White Sands method or the gradient refraction method, which allows accurate analysis of very low elevation 
angle and long-range data. Refraction properties are extrapolated from surface conditions, or a measured profile 
may be input. Velocity is determined by differentiating position. Accelerations are determined by differentiating 
velocity. This paper describes the algorithms used, gives the operational details, and discusses the limitations and 
errors of the program. Appendixes A through E contain the derivations for these algorithms. These derivations 
include an improvement in speed to the exact solution for geodetic altitude, an improved algorithm over earlier 
versions for determining scale height, a truncation algorithm for speeding up the gradient refraction method, and a 
refinement of the coefficients used in the White Sands method for Edwards AFB, California. Appendix G contains 
the nomenclature. 

INTRODUCTION 

One way to determine the trajectory of a flight vehicle is to analyze the ground-based, radar-tracking data. The 
radar measures range to the vehicle, azimuth of the vehicle from true north, and elevation angle of the vehicle 
above the local horizon. These measurements need to be filtered and corrected for atmospheric refraction. Then, 
geometric principles can be applied to convert these data into such familiar forms as latitude, longitude, and 
altitude. Next, derivatives can be taken to determine velocity and acceleration. These quantities are easily 
calculated in real time at the radar site, and the results are sufficiently accurate for many applications. To achieve 
these calculations in real time, however, certain assumptions are made regarding the structure of the atmosphere. 
These assumptions may introduce errors in the refraction corrections. Often for the high accuracies required for 
flight research, analyzing the atmospheric parameters after the flight and then using the results for refraction 
corrections is necessary. The process of analyzing the atmosphere can be quite involved 1 and will not be available 
in real time in the foreseeable future. Using atmospheric data gathered and analyzed before radar-tracking time 
suffers from temporal variations of the atmosphere. 1 An added benefit of analyzing the radar data after the flight 
using an atmospheric analysis is that such airdata parameters as Mach number, true airspeed, and pressure altitude 
may be accurately determined. 

A postflight FORTRAN program called “radar” reads and analyzes ground-based radar data from any radar site. 
This program provides Earth relative position, velocity, and acceleration parameters. Airdata parameters are also 
provided if atmospheric characteristics are input. This program reads data from the NASA Dryden Flight Research 
Center (NASA Dryden), Edwards, California, Flight Data Access System (FDAS) ; the encoded 9-track radar tape 
from the Army-Navy/Fixed Position System (AN/FPS-16) radars at Edwards AFB, California; or binary range, 
azimuth, and elevation angle data from any other radar. Cartesian position with respect to the center of the Earth 
can also be used as input. This position may be from an inertial navigation system or the Global Positioning 
System (GPS). Output from this program is in NASA Dryden compressed format. 2 As an option, the output can 
also be written in binary format. Program options include spike removal, data filtering, and atmospheric refraction 
corrections. Atmospheric refraction can be corrected by using either the quick White Sands method 3 , which is 
inaccurate at low elevation angles, or the accurate, but computationally slow, gradient refraction method. 4 

•Maine, Richard E., User’s Manual for GetFdas, Version 0.72, Apr. 30, 1993, NASA Dryden working paper. 
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Refraction properties are either extrapolated from surface conditions or determined from a user-supplied table of 
refraction as a function of altitude. The program models the Earth as an ellipsoid. Velocity is determined by 
differentiating position, and accelerations are determined by differentiating velocity. 

This paper describes the algorithms, operational details, limitations, and errors for version 6.2 of “radar.” 
Although earlier versions have existed for decades, this paper is the first time the program has been documented. 
Appendixes A through E contain the derivations for these algorithms. These derivatives include an improvement in 
speed to the exact solution for geometric ellipsoid altitude, an improved algorithm over earlier versions for 
determining scale height, a truncation algorithm for speeding up the gradient refraction method, and a refinement 
of coefficients used in the White Sands method for Edwards AFB, California. Appendixes A through E are 
universally applicable. Appendix F is specific to the program used at NASA Dryden. The nomenclature is given in 
appendix G. 

METHOD DESCRIPTIONS 

The methods for analyzing radar data are discussed next. These methods include geodetics, refraction 
corrections, spike removal and filtering, velocity and acceleration determination, and Earth relative and airdata 
parameters. 

Geodetics 

The basic information provided by a ground-based radar is time-referenced range, azimuth, and elevation angle 
to the vehicle. After various corrections are applied to these three quantities, the position of the vehicle with respect 
to the radar site can be calculated. Location of the radar site with respect to the center of the Earth can also be 
rftlrnlawi Adding these two vector positions yields the position of the vehicle with respect to the center of the 
Earth. The equations used to determine position are derived in appendix A. 

The Earth is modeled as an ellipse of revolution, otherwise known as an ellipsoid. Table 1 lists the semimajor and 
semiminor axes of this ellipse, a and b, in several systems. The World Geodetic Survey (WGS) 84 is used by most 
radar systems in the United States and is included in this table. Determining the altitude and latitude of the vehicle 
about an ellipsoid requires somewhat complex calculations, and these calculations are shown in appendix B. An 
im provement to the exact solution for determining altitude, which results in reduced computations, is also 
presented. 

Note that the algorithms in appendix B are used to calculate ellipsoid altitude which is different from geoid 
altitude. Geoid altitude is the altitude above mean sea level (m.s.l.). Because of mass irregularities of the Earth, the 
geoid is a highly irregular surface. For the continental United States, the geoid separation (ellipsoid altitude minus 
geoid altitude) is a negative number and is approximately -100 ft for radar 34 at Edwards AFB, California, using 
the WGS 84 system. Because most users desire geoid altitude, the ellipsoid altitude is biased by the radar site geoid 
separation to give geoid altitude. Another bias to altitude can be input to approximate the geoid separation change 
for flights great distances from the radar site. This method works well when the vehicle travels over an area where 


Table 1. Ellipsoid Earth models. 


Model 

Semimajor axis a, ft 

Semiminor axis b, ft 

al(a-b) 

WGS 84 

20925604.47 

20855444.88 

298.25722043 

WGS 72 

20925639.76 

20855480.71 

298.26002261 

“radar”, version 4.0 

20925832.00 

20854892.00 

294.97930764 


2 


geoid separation is relatively constant but is less effective for trajectories over long distances. A future 
improvement to this program would be an analytical calculation or database lookup of geoid separation as a 
function of latitude and longitude to adjust ellipsoid altitude to altitude above mean sea level. 

Refraction Corrections 

A radar unit measures the time a pulse of electromagnetic energy takes to travel from the radar antenna to the 
vehicle and return to its originating location. The initial assumption is that the pulse travels at the speed of light in 
a vacuum, so the range to the vehicle is easily calculated. The angle at which the antenna is pointed above the local 
horizontal is the measured elevation angle. The speed of light through the atmosphere is not the same as it is 
through a vacuum because it is affected by pressure, temperature, and humidity. Because these three quantities vary 
with altitude, the speed of light varies with altitude. The variation of the speed of light with altitude also causes the 
beam of light to bend. For example, consider two parallel beams of light in the atmosphere that are nearly 
horizontal with the Earth and at slightly different altitudes. The speed of light through the atmosphere is c / TJ, 
where c Q is the speed of light in a vacuum, and T| is the index of refraction. Because T| generally decreases with 
altitude, and c 0 is a constant, the upper beam travels slightly farther than the lower beam in the same amount of 
time. In this manner, the wave front has bent downward. This bending effect occurs for each incremental segment 
of the radar beam. 

Figure 1 shows the effect of a nonhomogeneous atmosphere on a radar beam. The radar beam follows a curved 
path, and most of the curvature occurs near the ground. The true straight line range to the vehicle is less than the 
measured range along the curved radar beam path. In addition, the true elevation angle is less than the measured 
elevation angle. This effect of atmospheric refraction is the greatest source of error and also the most difficult to 
correct. 

To correct for refraction, the properties of the atmosphere as a function of altitude must be determined. An easy 
method is to measure the properties at the radar site and then to extrapolate values at increased altitudes. Because 
of the extrapolation, this method is the least accurate. A more nearly accurate method is to measure the properties 


Spheres of 
constant 
refraction 
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/ vehicle 
/ location 


/ 


Apparent 

range 



True vehicle 
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True range 


Figure 1 . Effect of a nonhomogeneous atmosphere on a radar beam. 
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using a single weather balloon. The most nearly accurate method is to use all available weather data after the flight 
to generate a model of the atmosphere for the time of flight of the vehicle. 1 This method takes considerable effort. 
The first two methods may be used in real time, but the third method must be used after the flight. When the highest 
degree of accuracy is desired, the third method is used. 

Figure 2 shows a postflight-analyzed refractivity profile along with a profile extrapolated from surface 
conditions. The differences between the two models are significant for all but the shortest radar ranges. 

This computer program provides two methods of correcting for atmospheric refraction errors. The first is called 
the gradient refraction method, 4 which analyzes the radar beam one short segment at a time to determine the 
incremental bending between segments. This method is the most nearly accurate. On the other hand, because many 
small segments are analyzed at each time point, the method is computationally slow. The second method is called 
the White Sands method because it was developed at White Sands Missile Range, New Mexico. 3 This method uses 
an empirical fit to exact refraction corrections, such as gradient refraction results, at a given radar site. As a result, 
the coefficients used are geographically specific. The White Sands method is a function of radar site atmospheric 
parameters, so the structure of the atmosphere above the radar site is assumed. This assumption contributes the 
greatest error to the method. As an advantage, this method is computationally fast. 

Appendix C contains the derivations for index of refraction, gradient refraction, and White Sands methods and 
three improvements to the algorithms used in previous versions of “radar”. The first improvement deals with the 
algorithm for computing scale height. In turn, scale height is used to extrapolate index of refraction above the radar 
site. The algorithm used in years past induces as much as a 10-percent error in the atmospheric refraction 
corrections. 5 Appendix C also contains an alternate method that is substantially more nearly accurate. 5 

The second improvement is a new truncation algorithm for the gradient refraction method. This algorithm 
substantially reduces computation time required yet retains the accuracy of the method. Figure 3 shows examples 
of the percentage savings realized by this truncation algorithm as a function of elevation angle for two ranges. 

The third improvement deals with the White Sands method for Edwards AFB, California. The empirical 
constants used in years past for the White Sands method at NASA Dryden were inaccurate. New constants have 
been derived and are presented in appendix C. Figure 4 shows the error in vehicle position caused by elevation 

100 xIO 3 



Figure 2. Refractivity profiles from postflight-analyzed balloon data and extrapolated surface conditions, 
Vandenberg AFB, California, April 5, 1990. 
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Figure 3. Truncation algorithm savings for radar 34 N s = 338 and l s = 500 ft 


Error In position 
because of 
elevation 
correction 
error, 
ft 



Figure 4. White Sands method errors in position as a result of an elevation correction error, where zgeiod 
= 2662.6 ft, r m , = 600,000 ft, N s = 338, and l s - 500 ft. Here, the gradient refraction method is used as a truth 
model. 


correction error as a function of measured elevation angle. Here, the gradient refraction method was used as the 
truth model. Most flight test work with radar data at NASA Dryden has been conducted at elevation angles above 
10°. Refraction errors from using the White Sands method increase rapidly below this angle. Above 10°, using the 
new constants reduces the error in vehicle position by one-half as compared to using the old constants. 
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Spike Removal and Filtering 

Spike removal and filtering of the raw range, azimuth, and elevation angle data are important options of this 
program. The spike removal routine detects spikes as large jumps in the derivative of the data and removes them by 
using a hold-last-rate scheme to estimate the true value of the data. The user selects if spike removal is to be 
performed and sets the criterion for spike detection. 

Filtering of the raw data is performed with a discretized infinite impulse response (HR) filter 6 , which is a second- 
order, low-pass filter. The user selects the damping ratio and cutoff frequency of the filter. Appendix D provides a 
detailed description of the spike removal and filtering routines. 

Velocity and Acceleration Determination 

Once the position of the vehicle has been determined through geodetics, its velocity and acceleration can be 
calculated by taking derivatives. Such calculations are done by concatenating an open-loop zero with the IIR filter. 
The generation of noise associated with the taking of derivatives is suppressed by simultaneous low-pass filtering 
and differentiating. The user selects the cutoff frequencies for the velocity and acceleration filters. In addition, the 
user selects whether acceleration of gravity is subtracted from the vehicle acceleration. Acceleration of gravity is 
subtracted when radar acceleration is to be compared to acceleration from onboard accelerometers. Appendix D 
also gives the details of the velocity and acceleration calculations. When the cutoff frequency is selected to be zero, 
a second-order accurate, backwards difference differentiator is used instead of the DR filter differentiator. 7 

Earth Relative and Airdata Parameters 

The Earth relative parameters of total velocity, V, flightpath heading, ¥, and flightpath angle, y, are derived in 
appendix E. If atmospheric data as a function of altitude are input, then airdata parameters can be calculated. These 
airdata parameters consist of true airspeed, true Mach number, pressure altitude, Hp, ambient pressure, 
Ps^, and dynamic pressure, q, and are also derived in appendix E. The atmospheric data needed are pressure 
altitude, ambient temperature, windspeed, wind direction, and lateral pressure gradient magnitude and direction. 
Appendix F gives the form of the atmospheric data. 

PROGRAM USE 

Appendix F contains specific instructions for running the program on the SUN 600 computer (Sun 
Microsystems, Incorporated, Mountainview, California), including all input and output files and parameters. 

METHOD LIMITATIONS 

This section describes the limitations and expected errors of the “radar” program. The primary limitation 
involves memory. Errors of note include refraction correction, spike removal, filtering, and differentiation. 

Memory 

The refraction and atmospheric tables can each accommodate a maximum of 100 altitude breakpoints. The size 
of arrays for filtering and differentiation is 75,000, allowing at least 60 min of 20-sample/sec data to be analyzed in 
one run of the “radar” program. These array sizes can only be changed in the source code. 
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Refraction Correction Errors 


Because the true position of the vehicle is unknown, determining the errors due to refraction corrections is 
difficult. Several potential sources of errors in computing the refraction corrections exist. First, the value for 
refractivity at the radar site may have errors because of temperature and pressure measurement errors. Above the 
radar site, the relationship of refractivity with altitude may be assumed to be decaying exponentially. The result of 
such an assumption can be quite different from the true profile as shown by figure 2. If a profile of refractivity with 
altitude is determined by a weather balloon, balloon measurement errors will add to the refraction correction errors. 
Both of these profiles are assumed to be constant horizontally and in time. Again, such an assumption can be 
erroneous. 

The gradient refraction method is regarded as the most accurate method for calculating refraction corrections. 
Accuracy of this method depends on the length of the segment, and its optimum value depends on the roundoff and 
truncation errors of the computer. The accuracy of the refraction profile of the atmosphere also affects the accuracy 
of the gradient refraction method. Generally, one or two least significant bits (LSB) of angle error may remain after 
gradient refraction corrections are applied. The magnitude of the error depends on the quality of the refractivity 
profile. For the AN/FPS-16 radars at Edwards AFB, California, one LSB is approximately 0.0027°. There are 
17 bits in a 360° circle. The White Sands method is an approximation to the gradient refraction method. Figure 4 
shows typical differences between the two methods. 

Spike Removal, Filters, and Differentiation 

The spike removal routine needs to have spike-free data for the first few data points. In addition, depending on 
the value of a user-defined rejection criterion, some spikes may remain in the data, or some good data which are 
somewhat erratic may have been removed. The user should always run the program with the spike remover option 
turned off and compare those results to the results of running the program with the spike remover turned on. The 
low-pass HR filter tends to smooth out jumps in the data that may be real, such as an acceleration jump of an air- 
launched vehicle. Even though the time lag induced by the filtering is removed by time shifting the filtered 
parameters, some phase shifting of these data remains at the higher frequencies. These filters have start-up 
transients, especially when taking derivatives, and this fact should be considered when choosing data start times. 

Other Expected Errors 

The variability of the geoid separation with geography induces errors in the altitude above mean sea level. A bias 
may be applied to take into account the difference in geoid separation from the radar site to another location, but 
this approach does not totally address the problem. 

Table 2 lists typical errors in data from the NASA Dryden AN/FPS-16 radars. 8 Note that these estimates 
represent the errors that would be present if no corrections were made, except for the best possible manual 
alignments. Some of these errors are considerably less than those given for other installations because rigorous 
calibrations and maintenance are performed routinely at this installation. From this table and the discussion in the 
Refraction Correction Errors subsection, mislevel, solar heating, and refraction correction errors are the largest 
errors in the radar data. Mislevel readings are taken periodically and are kept within specifications, currently 10 arc 
sec. This program allows for finer correction to mislevel as well as makes refraction corrections, but other sources 
of errors are currently neglected. 
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Table 2. Typical estimated errors in the NASA Dryden AN/FPS-16 radar. 8 


Source 

Type 

Typical, mrad 

Typical, ft 

LSB values 

Thermal, range 

Noise 

— 

2.6 

0.41 

Thermal, angle 

Noise 

0.02 

— 

0.40 

R-f axis shift 

Bias 

— 

~0 

— 

Droop, el 

Bias 

0.03 

— 

0.61 

Orthogonality 

Bias 

0.02 

— 

0.41 

Mislevel 

Bias 

0.05 

— 

1.00 

LSB precision 

Noise 

0.03 

3.2 

0.50 

Solar heating 

Bias 

0.05 

— 

1.00 

Wind force 

Bias 

0.01 

— 

0.20 

Antenna unbalance 

Bias 

— 

-0 

— 

Servo unbalance 

Bias 

0.01 

— 

0.20 

Dynamic lag* 

Bias 

0.01 

— 

0.20 

Glint* 

Noise 

0.00 

— 

0.00 

Vertical deflection 

Bias 

0.02 

— 

0.40 

Earth model 

Bias 

0.01 

— 

0.20 


*Target at 150n. mi. 


CONCLUDING REMARKS 

A postflight FORTRAN program called “radar” that reads and analyzes ground-based radar data from any radar 
site has been presented. This program provides Earth relative position, velocity, and acceleration as well as airdata 
parameters if atmospheric characteristics are input. A general description of methods used, program use, input, 
output, and method limitations has been given. Detailed derivations of algorithms are given in the appendixes. 

In addition to documenting algorithms that have been used in earlier versions of this program, this report presents 
several new techniques and refinements. These techniques and refinements include an improvement in speed to the 
exact solution for geodetic altitude, an improved algorithm for determining scale height, a truncation algorithm for 
speeding up the gradient refraction method, and a refinement of coefficients used in the White Sands method for 
Edwards AFB, California. 

Dryden Flight Research Center 

National Aeronautics and Space Administration 

Edwards, California, May 28, 1993 
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APPENDIX A 


GEODETICS DERIVATION 


The Earth can be described as an ellipse of revolution with semimajor axis, a, semiminor axis, b, and axis of 
revolution through the North and South Poles (fig. A- 1 ). Considering a slice of the Earth that contains the axis of 
revolution (fig. A-2), a circle of radius a can be constructed. 


z 




Figure A-2. Geodetic geometry for a point on the ellipsoid. 
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Three angles are associated with any given point on the surface of the Earth: X c , X, and p. The geocentric 
latitude, X c , is the angle between (u, v), the center of the Earth, and the equitorial plane. The geodetic latitude, X, 
is the latitude shown on maps and is the angle that the local vertical line through (u, v) makes with the equatorial 
plane. Figure A-2 shows the reduced latitude, p. 9 From this figure, 


u = acos(p) and v = frsin(P) 


(A-l) 


To find the relationship between these three angles, the following equation of an ellipse is used: 



Differentiating implicitly and rearranging for the slope of the ellipse gives 


dv 

du 


,2 

-b u 


2 

a v 


The local vertical to the ellipse has a slope of 


= ^ = tan ^) 


(A-3) 


(A-4) 


In addition from figure A-2, 


tan(X ) = - 
L u 


(A-5) 


As a result, 


tan(X) = tan (A. ) 

b 2 


(A-6) 


Differentiating equation (A-l) gives 


du 

dv 


= f tan(P) 


(A-7) 
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Thus, 


a /jV 

tan (A.) = -tan(P) = z tan(A, ) 
b \b) 


(A-8) 


Having expressions for the sine and cosine of reduced latitude in terms of geodetic latitude will be useful later in 
this dicussion. Using equation (A-8) and the definition of the tangent gives 


tan(p) = 


-sin (A,) 
a _ 5 

cos (A,) C 


(A-9) 


Solving for the sine and the cosine, while noting that the eccentricity, e, of the ellipse is by definition 


e 2 =l-f b - 


(A-10) 


which gives 


Similarly, 


sin(P) = 


-sin (A,) 
a 


V^ + c 2 


sin 2 (A,) + cos 2 (A,) 


-sin (A.) 
a 


-sin (A.) 
a 


1(1 - e 2 ) sin 2 (A.) + cos 2 (A,) -e 2 sin 2 (A,) 


COS(P) = 


J^C 2 


cos (A.) 

^1 - e 2 sin 2 (A.) 


(A-ll) 


(A- 12) 


The position of a radar site is generally given in terms of geodetic latitude, A.^, longitude, 0 5 , and height above the 
ellipsoid, h s , where the s subscript denotes “of the radar site” From equation (A-l) and adding the increment for 
height, the position of the radar site in terms of u and v is 


u s = acos(p ) + h cos(A. ) 


(A- 13) 


Vj = fcsin(P ) + h sin(A, ) 


(A- 14) 


it 



Getting the position of the radar site in terms of Cartesian geocentric coordinates, xc, yc, and zc is desired. The x 
axis lies in the equatorial plane and points towards the prime meridian (0° longitude), the z axis points towards the 
North Pole, and the y axis completes the right-handed system as shown in figure A- 1 . Because 


xc = «cos(0) and yc = usin(0) 


(A- 15) 


using equations (A-l 1) and (A- 12), the Cartesian geocentric coordinates for the radar site are shown to be 


xc, = 


a 


+ h r 


yjl-e 2 sin 2 (X,) ) 


cos(A, )cos(0 ) 


(A- 16) 


yc* = 


a 


+ h„ 


%/l -e 2 sin 2 (X,,) ) 


cos (K s ) sin (0,) 


(A- 17) 


zc, = 


<*( 1-0 


■Jl-e 2 sin 2 (>t,) 


+ h . 


sin(X,) 


(A- 18) 


Now that the location of the radar site is known, the position of the vehicle with respect to the radar site can be 
added to this location to obtain the total vehicle position vector. 

Consider a right-handed coordinate system centered at the radar antenna that points locally north, east, and down, 
X[, yj, and Z[ . Assuming that the range, azimuth, and elevation have been corrected for various errors, the position 
of the vehicle is 


x t = r r cos(az)cos(e/ r ) 

y/ = r r sin(az)cos(el r ) (A-19) 


z, = -r r sin(el r ) 
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These components are transformed into the geocentric coordinates by rotating them through latitude and 
longitude and adding the radar site position to get 


cos(0,) -sin(0,) 0 
sin(0 J ) cos(0 J ) 0 
0 0 1 


cos(- 180° -X,) 0 

0 1 
-sin(- 180° - X s ) 0 


sin(- 180° -X,) 


x l 


* C 5 

0 


yi 

+ 


cos(- 180° -X s ) 


h 





-sin(X,)cos(0 5 ) -sin(0,) -cos(X,)cos(0,) 

x i 


xc. 


XC 

= 

-sin (X,) sin (0,) cos(0,) -cos(X 5 )sin(0,) 

yi 

+ 


= 

yc 


cos(X,) 0 -sin(X,) 

Z/ 




zc 


(A-20) 


Now that the geocentric coordinates of the vehicle are known, the altitude of the vehicle above the ellipsoid can 
be determined (appendix B). Part of this process involves finding the point on the Earth directly below the vehicle, 
known as the piercing point, x$, y^, z®, (fig. A-l). Because equations (A-l) through (A-12) are only valid for a 
point on the surface of the ellipse, the piercing point must be used to determine vehicle latitude. Using equation 
(A-5), the geocentric latitude is 


X c = tan 


-l 




/*© + y e J 


Geocentric latitude can be converted to geodetic latitude using equation (A-8) 


X = 


tan 



The longitude of the vehicle is determined by 


0 = 



(A-21) 


(A- 22) 


(A-23) 


Knowing the position of the piercing point relative to the radar site is often useful, such as when the ground track of 
the vehicle is desired. A difficulty arises on a nonflat Earth because traveling a certain distance north and then east 
will result in a different point than traveling first east than north, especially near the poles. For this reason, the 
distance from the radar site north to the latitude of the piercing point, xr, and the distance east from the radar site to 
the longitude of the piercing point, yr, is computed. 

The distance xr can be determined by taking the arc length along the ellipse of the Earth, that is, along the radar 
site longitude, from the radar site to the piercing point latitude and substituting equations (A-l) and (A- 10) to get 
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xr 


= J J(lp ) 2 + ("Spl dP = / </« 2 sin 2 (P) + » 2 cos 2 (P) dp 

(A-24) 

= a J Jl-e 2 cos 2 ((3) dp 

S 

which is an elliptical integral and must be calculated numerically. The “radar” program completes such calculations 
using Simpson’s Rule. A linear interpolation is used to select the values for a given p, and the difference in values 
from the radar site and piercing point is xr. 

The distance yr is the length of a circular arc from the vehicle longitude to the radar site longitude along the radar 
site latitude. From equations (A-l) and (A-12), 


yr = u s (e-Q s )-r£rz = acos^Xe-e,) 71 


acos(A,_) 


180 c 


180 c 


■T- 


2 . 2,* v 

e sin (A ) 




(A-25) 


Another geodetic quantity required for appendix C is the local radius of curvature of the Earth at the radar site. 
Re. From the definition of radius of curvature. 


Re = 


Ti f 

L 1+ U 


f diA 2 ! 1 - 5 


J J 


d 2 u 


dv 


(A-26) 


Taking the derivative of equation (A-7) gives 

d (-^ tan (P)) 


d 2 « 

dv 2 


dp _ 


-a 


dp dv 


b 2 cos 3 (P) 


(A-27) 


Substituting equations (A-7) and (A-27) into equation (A-26) gives 

1.5 


{1 +[|tan(P)] 


Re = 


[ a 2 sin 2 (p) + f > 2 cos 2 (P )] 15 ^> 2 cos 3 (P) 


a 


b 2 cos 3 (P) 


b 2 cos 2 (P) 


a 


(A-28) 


{a 2 sin 2 (P) + fc 2 cos 2 (P )} 1 5 

ab 
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From the definition of sin(P) and cos(P) in equations (A-l 1) and (A- 12), and noting that the radius of curvature at 
the radar site is required. That is, 


Re 


{ b 2 sin 2 (X 5 ) + b 2 cos 2 (X s ) 

1 l-e 2 sin 2 (^) 

ab 



l-e 


sin (X.) 



(A-29) 
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APPENDIX B 


ALTITUDE DERIVATION 


An exact method for determining altitude above an ellipsoid is presented by Hedgley 10 , and is outlined here with 
an improvement. In the original method, a fourth-order polynomial is solved. Each of the four roots must be used to 
determine four altitudes, and the lowest altitude is the correct answer. It is proven below that the minimum root is 
always the correct root. Because radar data are typically measured at 20 samples/sec and altitude is computed at 
every time point, this approach provides a substantial savings of computation time. 

The distance from the vehicle to a point on the ellipsoid is given by 

d = J(x - .xc) 2 + (y - yc) 2 + (z - zc) 2 (B-l) 


where (x, y, z ) is a point on the ellipsoid, and (xc, yc, zc) is the vehicle point. Figure A-l shows this geometry. The 
minimum distance is the altitude and is determined using the Lagrange multiplier method 11 where 


J(x, y, z) = (x - xc) 2 + (y - yc) 2 + (z - zc) 2 - a 


( 2 2 2 \ 

fL + 2L + 5__i 

u 2 a 2 b 2 ) 


Taking partial derivatives and equating them to zero gives 


dJ , cl2x n 

5 - = 2(x-xc) — = 0 

dx a 2 


or 


x = 


xc a 
2 

a - a 


Y = 2(y-yc)-^ = 0 

dy n 2 


or 


v = y^±- 
y 2 
a -a 


dJ , a2z ~ 

3 - = 2 (z-zc) = 0 

dz u 1 


or 


z = 


zc_b 
b 2 - a 


<LL = _£?_>L_z! + i = o 
« 2 « 2 b 2 


(B-2) 


(B-3) 

(B-4) 

(B-5) 

(B-6) 


Substituting equations (B-3) through (B-5) into equation (B-6) to eliminate x , y, z and collecting terms in a yields 


rj_v_ 2 ri + i 


r 2,2 

2 

2 

2l 

-la 3 + 

4 + 2- + L- 
,2 2 
L b a 

xc 

yc 


1 2,2j a 2 L 2 + , : 

L a b J L a b 

2 ja + 

b 2 

Ics 

1 

G 1 

ro| 

i 


(B-7) 


„r 2 2 2 2 ,2- r 2, 2 2,2 2,2 2 2, . 

+ 2 [xc +yc +zc -a -b ]a + [a b -xc b -yc b -zc a ] = 0 
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2 

A typographical error exists in the coefficient for a in Hedgley's paper. 10 Parameters a, b, xc, yc, and zc are all 
scaled to a to avoid using large numbers, then rescaled back once the altitude and piercing point have been 
determined. The original method uses each of the four possible real roots in turn to calculate altitude, and the four 
altitudes are compared to find the lowest altitude, which is the correct answer. 

The improvement reported here shows that the minimum real root gives the correct altitude. Taking equa- 
tion (B-3), noting that the point on the ellipse is the piercing point, and rearranging for a gives 


a = 


a 2 (x 9 -xc) 




(B-8) 


The vehicle and the point below it must be in the same hemisphere, so xc and x& have the same sign. For positive 
altitude, the magnitude of xc is greater than the magnitude of x$. As a result, a must be negative for positive 
altitude. For negative altitude, the magnitude of x© is greater than the magnitude of xc , so a must be positive. For 
zero altitude, a must be zero. Similar arguments hold using equations (B-4) and (B-5). 

Determining if any negative roots of equation (B-7), exist will be helpful. By Descartes' rule of signs 12 , the 
number of negative roots of /(a) = 0 is equal to the number of positive roots of /(-a) = 0 , which is equal 
to the number of sign changes of the coefficients of / (-a) or that number reduced by a positive even number. In 
equation (B-7), normalize all variables by a, so a' equals one, b' is slightly less than one, and the sum of the squares 
of position components are nearly equal to the square of the distance from the center of the Earth, R a ’ 2 . With these 
assumptions, using -a in place of a to look at the signs of the coefficients gives 


[ + ]a 4 + [ + ]a 3 + [-6- - R 0 2 ]a 2 + 2[ - 2 - R 0 ' 2 ]a + [ - 1 - - R 0 ' 2 } = 0 (B-9) 


For the case of -1 > ~R 0 ' 2 (which is negative altitude), the signs of equation (B-9) are all positive, so no sign 
changes. For this reason, no negative roots exist for negative altitude. In all other cases, only one sign will change, 
so only one negative root exists for positive altitude. 

Because xc and x© have the same sign, a must be less than a from equation (B-3). From equation (B-5), a must 
be less than b . Because for Earth, b is less than a, a must be less than b . Now, substituting equations (B-3), 
(B-4), and (B-5) into equation (B-l) gives 


d 


2 


(xc + yc) 


2 


a 

2 2 

r « i 

2 

a -a 

+ zc 

i — 
1 


(B-10) 


2 

Figure B-l shows a graph of d as a function of a. Note that this figure shows one negative root for positive 

altitude and no negative roots for negative altitude. To prove that the shape of these curves is correct, the derivative 
with respect to a is taken. That is. 


did 1 ) 

da 


2 

2 (xc + yc) 


2 1 

r 


,2 ^ 

r 

a 

a 

+ 2zc 2 

b 

a 

2 2 
{(d -a) J 

2 

L a - a_ 

2 2 
{(b 2 -a) J 

,2 

b -a. 


(B-l 1) 
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Figure B-l . Equation (B-10) roots as a function of Lagrange multiplier. 


Only the terms in square brackets affect the sign of the derivative. For a positive, and noting that a is less than 
, the derivative is positive. The larger a is, the bigger d 2 becomes. For a negative, the derivative is negative. 
Again, the larger the magnitude of a is, the bigger d 2 becomes. This fact is visualized by the concave upward shape 
of figure B-l. Because the root that gives the smallest value of d 2 is desired, the minimum positive real root or the 
maximum negative real root yields the true altitude. 

For negative altitude, all the roots are positive. If the roots are all positive, then the minimum root yields the tme 
altitude. For positive altitude, only one negative root exists, and the correct root must be negative. As a result, the 
minimum real root yields the true altitude; therefore in all cases, the minimum real root of equation (B-7) gives the 
correct value for altitude. 

Now, determining the roots of a fourth-order polynomial is required. The solution to a quartic is given by 
Dickson 12 and is based on the work of Ferrari (1522-1565). Given quartic 


* f 4 + V * 3 + c q x q + d q x q + e q = 0 


rearrange to get 



+ b„ x. 
q q 


- -c„x„ 

q q 


- d„x„ - e„ 

q q q 


(B-12) 


(B- 13) 
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Complete the square on the left-hand side 


b„ \2 (b 


c„ x„ -d„x„-e„ 
9 9 q q q 


Adding 


( 2 °q \ y c 

x. + -Jjc v, + — — 

y q 2 9) J c 4 


to both, sides gives 


2 b y c 
x„ + -£x. + — 
q 2 9 2 


, 2 \ , 2 

b q 2 ( b q y c y c 

—L — c„ + V. x„ + d n x„ + — — e„ 

4 q •'c q y 2 9 ) 9 4 9 


and rearranging to get the resolvent cubic equation 


v 3 - c n y} + ( b„d n - 40 v - (b 2 e„ + 4 c n e. - d 2 ) = 0 
J c qJ c v q q q'^c v q q q q q ' 


y c + b c y c + C c y c + ^ c 


Now find any real root y c of the resolvent cubic equation by setting 


•Vc = *e"T 


giving 


z c +pz c + q = 0 


which is the reduced cubic equation with 


b c b c c c 2b c 

p = c, and a — d + 

r C 7 v "c 'J '77 
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The discriminant of the reduced cubic and of the resolvent cubic is 


A = -4p 3 -21q 2 


(B-21) 


If A<0 , then one root is real and two are complex. If A = 0 , then all the roots are real, and two or more are 
repeated. The real solution for A < 0 is given by 


If A > 0 , then three distinct real roots exist, and the following trigonometric solution is used: 

cos(3Z) = 4 cos 3 (Z )-3cos(Z) 

Replacing Z by Z + 120° and Z + 240° , in turn, gives 

cos(3Z+ 360°) = cos(3Z) = 4 cos 3 (Z+ 120°)-3cos(Z+ 120°) 
cos(3Z+ 720°) = cos(3Z) = 4 cos 3 (Z+240°)-3cos(Z+240°) 

Thus, cos(Z ) , cos(Z + 120°) , and cos(Z + 240°) are the three roots of the equation 

4/ c 3 -3 t c = cos(3Z ) 

Hence, 

. 3 _ 3 _ cos(3Z) 0 
*c 4 *c 4 u 

To solve the reduced cubic equation, take z c = st c . The result is 


(B-22) 


(B-23) 


(B-24) 


(B-25) 


(B-26) 


s s 


(B-27) 
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with 


s = 



cos(3Z> = f(|)' 5 


As a result. 


z 


c 



(B-28) 

(B-29) 


(B-30) 


Now that z c is known, a root of the resolvent cubic equation can be found. That is. 




(B-31) 


Returning to the quartic equation, the right-hand side of equation (B-16) is the square of a linear function. For 
example, mx q + n . Thus, 



9 

a2x +b2x +c2 = 
<1 ^ 


2 

(mx q + n) 


(B-32) 


In addition. 




mx q + n 


or 




= -mx q + n 


(B-33) 
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If a 2 = 0 , then the polynomial is a perfect square. The four roots to the quartic are the four roots of the following 
quadratics: 



b q Je 
—x + — 
2 * 2 


-Je 2 = 0 


(B-34) 



+ 


y_c 

2 


+ Jc2 - 0 


If a2*0, then the polynomial is not a perfect square. The four roots to the quartic are the four roots of the 
following quadratics: 


(B -35) 

V + (f-'”K + T + " = 0 


where 


m = Jal 


n = 


bl 

2ja2 


(B-36) 


When solving the two quadratic equations, the positive of the radical may be neglected because only the minimum 
roots are of interest. As a result, there will be one or two real roots from which to choose the minimum root. 

Once the correct a is determined (the minimum a), the piercing point can be found by using equations (B-3), 
(B-4), and (B-5). The geocentric latitude, geodetic latitude, and longitude of the vehicle can be determined from 
equations (A-21), (A-22), and (A-23). Rearranging equations of the form of equations (A-16), (A-17), and (A-18) 
allows solving for the altitude of the vehicle by using one of the following three equations: 


h' = 


h' = 


xc 


cosWcOS(e) 7l-« 2 sin 2 W 


yc 


cos (X) sin (0) f 2 . 2.* v 

VI -e sin (A.) 


(B-37) 

(B-38) 


W = -T 


zc 


1-0 


sin (A.) f 2 . 2 ... 

Vl -e sm (A.) 


(B-39) 
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where 


h = h’ a (B-40) 

The proper equation is chosen to avoid division by values close to zero in the first term. These equations determine 
altitude with a higher computational accuracy than can be obtained by using equation (B-l). Lastly, because all 
parameters were scaled to a to avoid using large numbers, the altitude is scaled back by the value a. 
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APPENDIX C 


REFRACTION CORRECTION DERIVATION 


To calculate the amount of bending in the radar beam, the index of refraction of the atmosphere must be 
determined. The index of refraction is T| = 1.0 in a vacuum, and T| > 1.0 in the atmosphere. This index decreases 
with altitude in most cases. When dealing with air, it is useful to deal with refractivity, N, where 

N=( ri-l)10 6 (C-l) 

The refractivity ranges from 0 in space to the order of 300 at the ground; therefore, t| ranges from 1.000000 to 
- 1.000300. The index of refraction at the radar site can be determined from site pressure and temperature 
measurements. Above this site, the index of refraction can be extrapolated or measured using weather balloons. 

Refractivity at the radar site, N s , is determined using the dry and wet bulb temperatures (°R), T s and Twet s , 
and the ambient pressure (in. Hg.), p s , at the site by 13 


N s = 


4730.3 341.36ev 4.1146 x 10 ev 
T T T 2 

1 s s 1 s 


(C-2) 


where 14 


b 

ev = Twet* 10 


b 

Twet,j 


•1 


[f + g(Twet s - d )]p s (T s - Twet s ) 


es = 




Table C-l gives constants a through g. The relative humidity given by 

rh = — X 100 

es 


(C-3) 

(C-4) 


(C-5) 


Table C-l . Constants for equations (C- 3) and (C-4). 14 


Constant 

Twet s above freezing 

Twet s below freezing 

a 

-4.9283 

-0.32286 

b 

-5287.32 

-4869.38 

c 

23.2801 

10.0343 

d 

459.4 

459.4 

f 

3.595 x 10- 4 

3.595 x lfr 4 

g 

2.336 x 10- 7 

2.336 x 10- 7 
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is checked to ensure that these pressures and temperatures do not yield values above 100 percent. Equation (C-2) is 
valid to within 0.5 percent for temperatures from -58 to 104 °F, pressures from 5.91 to 32.48 in. Hg., ev from 0 to 
0.88 in. Hg., and radar frequencies up to 30,000 MHz. 13 

One way of determining refractivity at higher altitudes is to assume that it decays exponentially with geoid 
altitude. That is. 



where h s is the ellipsoid altitude of the radar site, gs s is the geoid separation at the radar site, ze is the ellipsoid 
altitude of interest, and H is the scale height. 

Scale height can be assumed by iterating a function of N s , using h s to determine function coefficients. 5 
The coefficients of this function were determined from a least squares fit of refraction correction data from nine 
tracking sites located throughout the world. The scheme uses an initial estimate of H = 7000m and then iterates 
the equation 



until H converges where the coefficients A, B, and C are given in table C-2. 5 The program defines convergence as 
within 1 ft. This method is superior to the one previously used at NASA Dryden. 4 - 15 The former method can cause 
errors in excess of 10 percent in radar refraction corrections for such semiarid areas as Edwards AFB, California. 5 

Another approach to determining refractivity at higher altitudes is to take a profile of refractivity as a function of 
altitude as measured by weather balloons. This approach is especially appropriate when the atmosphere has a high 
degree of nonexponentiality, such as when an inversion layer is present near the surface where most of the radar 
beam bending takes place. When interpolating between data points, altitude is interpolated linearly, and N is 
interpolated exponentially. The “radar” program extrapolates using the nearest two points for altitudes above and 
below the profile, so the end points should exhibit the same trend as the rest of the data. 

This model of atmosphere refractivity is static, but the “radar” program could be modified to accept a time 
history of refractivity to do dynamic atmospheric refraction corrections. To date, no attempt has been made to take 
into account horizontal gradients of refractivity of the atmosphere. This effect may become significant during 
atmospherically active days, but those days would hopefully see little flight activity. With these caveats, refractivity 
can be determined for any altitude; now, a determination of how it affects the radar beam can be made. 


Table C-2. Constants for equation (C-7). 


Constant 

h s < 1000 m 

1000 m<h s < 2500 m 

h s t 2500 m 

A, m 

17590.00 

18588.000 

21273.000 

B, m 

30.55 

40.814 

60.227 

C, m 

0.00 

1500.000 

3000.000 
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GRADIENT REFRACTION METHOD 


For the gradient refraction correction, the time the radar beam takes to travel from the radar to the target is 
divided into a number of equal parts, delt being the time interval, and ns being the number of segments. It is 
agcnnv-H that the radar beam travels in a straight line during each time interval and has a discrete angular change 
between adjoining segments. The first segment leaves the radar at the measured elevation angle (fig. 1). Now 
assume that each segment consists of three beams (fig. C-l). The time interval in which the three beams travel is the 
same. Because the index of refraction, and thereby the speed of light, changes with altitude, the top beam travels 
the farthest, and the lower beam travels the least distance. As a result, the wave front becomes increasingly vertical, 
and the width of the beam increases. 

Now, consider two adjoining segment triplets to see how the turning is calculated. The nth triplet is some distance 
from the radar site, and the local vertical there is tilted by the internal Earth angle, ia n , from the radar site vertical. 
The beams Us n , Ms n , and Ls n are at an elevation el n + ia n to the local horizontal. Beam Ms n starts at an 
altitude of h n . Next estimate the ellipsoid altitude of the midpoint of beam Ms n as 


c„ sin(e/_ + /fl_) 

7 o j j. v n n / 

ze n = H delt - 


" n„_i 


(C-8) 


where c Q is the speed of light in a vacuum, and T| n _ , is the index of refraction from the previous segment. Now, 
the refractivity at altitude ze n , N n , can be determined either through table lookup of balloon data or through equa- 
tion (C-6). 

The derivative of refractivity with respect to altitude, ^ , can be determined by differentiating Equation (C-6). 
That is, 



3937 ft 
1200 m 


(C-9) 
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where H is determined by equation (C-7) or by adjacent points in the balloon data. In the case of balloon data, H 
may differ for each segment. The index of refraction for the upper and lower beams is now 


N n, = N n + (t^) n K '" C0S(W » + “O (C-10) 

- g^)^„cos( e ;„ + la.) (C-ll) 

where w n is the distance between the midbeam and the upper or lower beam of the current segment. 

The lengths of the three beams are as follows: 


Ms n 


—delt 
— delt 
—delt 


(C- 12) 
(C- 13) 
(C- 14) 


The width of the next segment is 


w 


n + 1 



(C-15) 


The turning angle is 


6 


n 


sin 


( 


2w n+ 1 J 


(C-16) 


A plane tangent to the radar site can be defined, where D is the distance downrange of the radar site, and Z is the 
altitude above the tangent plane. The increment to tangent plane altitude and downrange because of this segment 
can be calculated in feet by 


D n + 1 = D n + Afs n cos(e/ n ) (C-17) 

Z n+ \ = Z n + Ms n sm(el n ) (C-18) 
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Note that the angle ia is not used in equations (C-17) and (C-18) because a radar site origin is being used instead of 
one local to the segment. The tangent plane position can be converted to altitude about a spherical Earth using the 
law of cosines and figure C-2. 


h n + 1 = A /(/?e + /t 5 + Z„ + 1 ) 2 + D n+1 -Re 


(C-19) 


To use equation (C-19), the radius of the Earth at the radar site. Re, is needed. Because the Earth is modeled as an 
ellipsoid, the local radius of curvature of the ellipse will be used. This radius is given by equation (A-29). The 
radius of curvature only needs to be computed once outside the refraction calculation loop because this radius only 
depends on radar site geodetic latitude and Earth characteristics. As an alternative to equation (C-19), the exact 
calculation for altitude given in appendix B may be used although this calculation greatly increases the amount of 
computation time required. 

Next, the internal Earth angle, ia, is calculated. For the spherical Earth, this angle is obtained by the definition of 
the sine 


ia 


n + 1 



(C-20) 


For the ellipsoid, the definition of the dot product is used. The two vectors are the site vertical vector, s, and the 
local vertical vector, v. Both vectors are normal to their respective surfaces of the Earth; therefore. 


ia 


n + 1 


-1 

COS 




V 


+ ^2 v W -H 2 + 5 3 V n + l 3 > 

*|K + l| > 


(C-21) 



Figure C-2. Gradient refraction method geometry for a spherical Earth. 
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Once, the parameters have been calculated for a single segment, the process is repeated for the subsequent seg- 
ments with equations (C-8) through (C-19) and either (C-20) or (C-21) until all the segments have been analyzed. 
Referring back to figure C-2, the corrected range and elevation are determined by 


r r = JdJ 

el r = tan -1 



(C-22) 

ft) 

(C-23) 


The accuracy of the gradient refraction method depends on the length of the time segment, delt, and its optimum 
value depends on the roundoff and truncation errors of the computer. For the SUN 600 computer at NASA Dryden, 
a value of delt that gives 1000-ft segments in a vacuum is generally used. Although this value may not be the opti- 
mum, it was selected because data from balloon profiles of reffactivity are received in 1000-ft intervals. 

Because the largest gradient of reffactivity with altitude occurs near the Earth, most of the bending of the radar 
beam occurs near the Earth. Above a certain altitude, the remaining turning in the radar beam may be insignificant. 
If the algorithm can be truncated above this point, a substantial computational savings is possible. An original 
method is presented now to determine when the gradient refraction algorithm may be truncated without introducing 
significant errors. 

The turning angle decreases with increasing altitude. This decrease is assumed to be at the same rate as 
reffactivity; therefore, the approximated turning angle is 


5 


n 


ft 


„ 3937/f 

1200m- 1 


(C-24) 


where A h n is an estimate of the altitude of the target above the current segment 


Ah n = (ns-n+ l)Ms n sin(el n + ia n ) 


(C-25) 


The difference between the free-space segment length and the segment length at a given altitude is likewise approx- 
imated as 


-A*. 

,, 3937/t 

. i onn m 

(c 0 delt-Ms n r = ( c 0 delt-Ms n )e 12 °° m (C-26) 
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The average turning that remains for each segment and the average segment length can be found by integrating 
with respect to altitude and dividing by the remaining altitude 


J 8 "dh \ S„e 


,, 3937/f 
L 1200m- 1 


8„ = 

r§ r+m 


_ K 


dh 3937/f g 
1200m " 


Aft, 


Aft, 


Aft, 


l-e 


„3937// 

L 1200m-* 


'•ns 

J j^c 0 c?e/f- (c 0 tfe/f-Mj rt )"Jrfft 


= 


Aft, 


= c 0 delt- 




z' 


Aft. 


-AA, 




l-e 


3937 /r 
1200m- 1 


(C-27) 


(C-28) 


The remaining turning is the average turning multiplied by the number of remaining segments 


s„.._ = 1 ) = 




Aft. 


-A/i. 


1/ 


l-e 


3937 /f 
1200m 


(C-29) 


The algorithm ends if the remaining turning is less than the minimum turning. Minimum turning is arbitrarily 
defined as 40 percent of the least significant bit of elevation. Because the digitization of elevation is 17 bits in a cir- 
cle, the minimum turning is 

360° 

8 m in = 0.4-221- (C-30) 

2 

If the algorithm is truncated early, the downrange distance from the radar site, D, and the altitude above the tangent 
plane, Z, need to be adjusted for the remaining segments. This adjustment is done by adding the components of the 
average segment length times the number of remaining segments. That is, 

D ns = D n + W n (ns-n+ l)cos(el n ) (C-31) 

z ns = Z n + M n (ns-n+ l)sin(e/„) 
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(C-32) 


When inversion layers and other atmospheric phenomena are present near the ground, the index of refraction does 
not decay exponentially. Lack of such exponential decay can invalidate equations (C-24) and (C-26). For this rea- 
son, the gradient refraction algorithm does not truncate below a certain critical altitude where these conditions 
might exist. The radar program defaults to a critical geoid altitude of 10,000 ft; however, the user may select a dif- 
ferent value. Figure 3 shows examples of the percentage of savings realized by this truncation algorithm as a func- 
tion of elevation angle for two ranges. 

WHITE SANDS METHOD 

The White Sands method for refraction correction was created out of a need to process radar data easily and in 
nearly real time. 3 This method uses an empirical fit to exact refraction corrections, such as results of the gradient 
refraction methods at a given radar site. For this reason, the coefficients used are geographically specific. Because 
only radar site atmospheric conditions are quickly available, the structure of the atmosphere above the radar site is 
assumed. This assumption contributes the greatest error to the method. The method was designed to provide 
accurate results for elevation angles from 1° to 90° and for ranges of 500 to 200,000 yd. 

A separate correction exists for elevation angle and range. Elevation angle correction will be discussed next. The 
refiractivity at the radar site is used to calculate the constant K le . 


y i n -6 6400 j. 

K " = 10 2 K 

(C-33) 

where there are 6400 army mils/2jt radians. Then, the measured downrange and vertical distance from the radar site 
in yards, D and Z respectively, are calculated by 

D = jcos (el m ) 

(C-34) 

Z = ysin (elj 

(C-35) 

These distances allow the elevation error in army mils to be calculated by 


K,,D 
Agl ~ K 2e + Z 

(C-36) 


where the constant K 2e will be discussed shortly. Lastly, the corrected elevation angle comes from 


el 


r 


el el 

m 6400 


(C-37) 
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Range is corrected in a similar manner by 


A r 


KlrD 

K 2r + Z 


so 


r_ = r -3Ar 

r m 


(C-38) 


(C-39) 


The constants K . 2e , K[ r , and K 2r are determined by a least squares fit of a large set of exact refraction correc- 
tions, over a range of values of and r m , for the desired radar site. These constants are a function of . The 
constant K 2e is given by 


*2e = 


K^DAel-^ZAel 2 


(C-40) 


The K lr and K lr are determined by the simultaneous linear equations (which are in error in reference 3) 


K lr ^D 2 - K 2r ^DAr = ^DZAr (C-41) 

K lr ^D Ar - K^Ar 1 = £ZAr 2 (C-42) 

which transform to 

(S DAr )fZ ZA '- 2 l - (I DZA ofI A '- 2 l 

K. r = i (C-43) 

d« A 0 2 -(I^ 2 )(Z Ar2 ) 

(Z D2 ¥Z ZAr2 l - (Z» ZA o(Z DAr ) 

K lr = 4 & (C-44) 

( Xd a o 2 -(I« 2 )(Z a ^ 

Table C-3 shows values of these four constants for radars at Edwards AFB, California, and White Sands Missile 
Range, New Mexico. 3 The values labeled Old Edwards have been used at NASA Dryden for several decades; 
however, no documentation for them exists. Values for K lr and K 2r did not exist for Edwards AFB, California, 
during that time. The values labeled New Edwards were computed using the gradient refraction algorithm with 
segment lengths of 500 ft; for ranges of 1,500, 3,000, 6,000, 15,000, 30,000, 60,000, 150,000, 300,000, and 
600,000 ft; and for elevation angles of 2°, 5°, 12°, 25°, and 70°. These ranges and elevation angles are also used to 
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determine the White Sands constants. Computations using the New Edwards constants have approximately one half 
the error that the Old Edwards constants yield for elevation angles above 10° (fig. 4). The potential exists for 
optimizing the constants for a given radar coverage by selecting certain combinations of range and elevation 
angles, but this effort is left for future work. 


Table C-3. Constants for equations (C-36) and (C-38). 


Old New Edwards White Sands 3 

Edwards zgeoid 5 = 2662.593 ft zgeoid s = 4000 ft 




K 2e , yd 

* 2 «. yd 

K lr , yd 

*2r> Yd 

K le , yd 

yd 

*2r’ yd 

220 

0.2241 

21410.6 

19915.6 

-3.627 

15789.2 

21790 

-3.023 

11980 

222 

0.2261 

21160.5 

19773.1 

-3.630 

15655.8 




224 

0.2281 

20916.6 

19630.4 

-3.633 

15522.3 

21200 

-3.025 

11820 

226 

0.2302 

20679.8 

19487.4 

-3.636 

15388.7 




228 

0.2322 

20449.3 

19344.2 

-3.638 

15255.1 

20650 

-3.029 

11680 

230 

0.2343 

20224.7 

19200.7 

-3.639 

15121.4 




232 

0.2363 

20005.9 

19056.9 

-3.640 

14987.7 

20130 

-3.033 

11540 

234 

0.2384 

19792.9 

18912.8 

-3.640 

14853.9 




236 

0.2404 

19585.2 

18768.4 

-3.640 

14720.0 

19610 

-3.037 

11400 

238 

0.2424 

19382.8 

18623.7 

-3.639 

14586.0 




240 

0.2445 

19185.4 

18478.7 

-3.638 

14451.9 

19110 

-3.041 

11270 

242 

0.2465 

18993.0 

18333.4 

-3.636 

14317.8 




244 

0.2485 

18804.8 

18187.8 

-3.634 

14183.5 

18650 

-3.046 

11140 

246 

0.2506 

18621.5 

18041.8 

-3.631 

14049.1 




248 

0.2526 

18442.3 

17895.5 

-3.627 

13914.7 

18250 

-3.051 

11020 

250 

0.2546 

18267.3 

17748.8 

-3.623 

13780.1 




252 

0.2567 

18096.4 

17601.7 

-3.618 

13645.3 

17900 

-3.055 

10890 

254 

0.2587 

17929.1 

17454.2 

-3.613 

13510.5 




256 

0.2608 

17765.9 

17306.3 

-3.607 

13375.4 

17550 

-3.059 

10760 

258 

0.2628 

17606.1 

17158.0 

-3.601 

13240.3 




260 

0.2648 

17449.7 

17009.3 

-3.594 

13104.9 

17200 

-3.064 

10640 

262 

0.2669 

17296.8 

16860.2 

-3.586 

12969.4 




264 

0.2689 

17147.1 

16710.5 

-3.578 

12833.8 

16870 

-3.069 

10520 

266 

0.2709 

17000.6 

16560.5 

-3.569 

12697.9 




268 

0.2730 

16857.2 

16409.9 

-3.560 

12561.8 

16550 

-3.074 

10400 

270 

0.2750 

16716.7 

16258.8 

-3.550 

12425.5 




272 

0.2771 

16579.0 

16107.1 

-3.539 

12289.0 

16250 

-3.079 

10280 

274 

0.2791 

16444.2 

15954.9 

-3.528 

12152.2 




276 

0.2811 

16312.0 

15802.2 

-3.517 

12015.2 

15970 

-3.085 

10170 
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Table C-3. Concluded. 




Old 

New Edwards 



White Sands 3 




Edwards 

zgeoid s 

= 2662.593 ft 

zgeoid s = 4000 ft 


278 

0.2832 

16182.2 

15648.8 

-3.504 

11877.9 




280 

0.2852 

16055.2 

15494.8 

-3.491 

11740.3 

15670 

-3.090 

10060 

282 

0.2872 

15930.4 

15340.2 

-3.478 

11602.4 




284 

0.2893 

15808.2 

15184.9 

-3.463 

11464.2 

15400 

-3.095 

9940 

286 

0.2913 

15688.1 

15028.9 

-3.449 

11325.7 




288 

0.2934 

15570.5 

14872.1 

-3.433 

11186.8 

15110 

-3.100 

9840 

290 

0.2954 

15454.9 

14714.6 

-3.417 

11047.5 




292 

0.2974 

15341.4 

14556.3 

-3.400 

10907.7 

14850 

-3.105 

9730 

294 

0.2995 

15229.9 

14397.2 

-3.382 

10767.6 




296 

0.3015 

15120.4 

14237.2 

-3.364 

10627.0 

14950 

-3.111 

9630 

298 

0.3035 

15012.8 

14076.3 

-3.345 

10485.9 




300 

0.3056 

14907.2 

13914.4 

-3.325 

10344.3 

14350 

-3.117 

9550 

302 

0.3076 

14803.3 

13751.4 

-3.305 

10202.1 




304 

0.3097 

14701.2 

13587.5 

-3.284 

10059.3 

14110 

-3.122 

9470 

306 

0.3117 

14600.9 

13422.4 

-3.262 

9915.9 




308 

0.3137 

14502.2 

13255.8 

-3.239 

9771.6 

13900 

-3.128 

9390 

310 

0.3158 

14405.1 

13089.1 

-3.216 

9627.5 




312 

0.3178 

14309.7 

12919.7 

-3.191 

9481.5 

13680 

-3.135 

9320 

314 

0.3198 

14215.7 

12749.5 

-3.166 

9335.1 




316 

0.3219 

14123.2 

12577.0 

-3.140 

9187.2 

13470 

-3.141 

9260 

318 

0.3239 

14032.4 

12403.5 

-3.113 

9038.7 




320 

0.3259 

13942.8 

12228.1 

-3.085 

8889.2 

13260 

-3.148 

9200 

322 

0.3280 

13854.7 

12051.0 

-3.056 

8738.5 




324 

0.3300 

13767.9 

11871.8 

-3.026 

8586.5 

13070 

-3.155 

9140 

326 

0.3321 

13682.5 

11689.8 

-2.995 

8432.6 




328 

0.3341 

13598.3 

11505.9 

-2.963 

8277.6 

12900 

-3.162 

9080 

330 

0.3361 

13515.4 

11319.6 

-2.930 

8121.0 




332 

0.3382 

13433.8 

11130.4 

-2.895 

7962.6 

12750 

-3.169 

9020 

334 

0.3402 

13353.2 

10938.2 

-2.859 

7802.1 




336 

0.3422 

13274.0 

10741.9 

-2.822 

7638.8 

12600 

-3.177 

8960 

338 

0.3443 

13195.8 

10542.4 

-2.783 

7473.4 




340 

0.3463 

13118.8 

10338.8 

-2.743 

7305.2 

12520 

-3.184 

8910 
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APPENDIX D 


SPIKE REMOVAL, FILTERING, AND DIFFERENTIAL DERIVATION 

The spike removal, filtering, and differentiation routines for the “radar” program are derived in this appendix. 
First, the spike detection and removal routines are described. Next, the transfer functions for the filtering and 
differentiation routines are presented. Then, the filters are discretized to give time recursive difference equations. 
Finally, the calculation of velocity and acceleration components is discussed. 

SPIKE DETECTION AND REMOVAL 

Automated spike detection and removal is an important utility of the program. Spike removal may be turned on 
or off. By default, such removal is not performed during the nominal code operation. The traditional difficulty with 
automated spike detection is in determining what constitutes a spike. Clearly, large data spikes can be detected and 
removed using such traditional statistical techniques as 3a. Unfortunately, the 3a technique will not reliably work 
for detecting subtle spikes which occur within the standard deviation of the data. 

The spike detection and removal routines implemented in this program overcome this difficulty by differentiating 
the suspected spike-corrupted signal. Differentiating the signal greatly amplifies data spikes and renders them 
clearly distinguishable from the input data stream. This effect is illustrated in figure D-l where a subtle spike in 
the input data stream is rendered extremely large in the differentiated signal. Clearly, if the original signal were 
used to perform the spike detection, the 3 a criterion would not have been violated. The spike would have remained 
undetected. 



Figure D-l. Effect of differentiation of input data spike. 
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Conversely, the differentiated signal clearly violates the 3a criterion at the data spike, and the data spike is easily 
detected. A sliding window, whose length is user defined, is used to accumulate sample mean and standard 
deviation statistics of the differentiated signal. The user also selects the number of standard deviations for the 
rejection criterion. When a data point is encountered whose magnitude deviation from the current sample mean 
exceeds the rejection criterion, then the point is rejected. The previous value of the derivative is used to extrapolate 
the original signal across the corrupted region. As such, the spike-editing routines perform a hold-last-rate 
interpolation. Mean and standard deviation statistics are calculated for a window twice: first using all derivative 
data, then again using derivative data that are within the criterion band. For spike detection, a first-order-accurate, 
backward-difference derivative is used. A single data spike will cause a large derivative at that point and another 
large derivative of opposite sign at the next point. As a result, two points will be removed. The filtered derivative 
described in the Differentiator Transfer Function subsection was not used. The effect of the one data spike would 
have been spread out over many data points because of the filtering. Hence, many data points would be removed. 

The spike detection routine demands that the first few data points be spike free, and the last one-half window of 
data points will be discarded. The spike remover should never be used blindly on new data. First, run the data 
without using the spike remover to determine if it is needed. If so, use various numbers of standard deviations as 
the detection criterion, and inspect the results. A high number of standard deviations may allow many spikes to go 
undetected, and a low number may remove valid data points. 

LOW-PASS FILTER TRANSFER FUNCTION 


The filtering and differentiation routines for the program are based on the use of an infinite impulse response 
(DR) filter to eliminate noise above a selected cutoff frequency. 6 In these routines, a second-order low-pass filter is 
used. The frequency domain transfer function of this filter is given by 


0(5) 

1(5) 



(D-l) 


where (d n is the natural resonance frequency, ^ is the damping ratio, and $ is the Laplace transform variable. For 
sinusoidal inputs, s = j (£> , the magnitude and phase angle of the filter transfer function may be written as follows: 


In addition, 


M(j to) = 





(D-2) 


(D-3) 


36 


Figure D-2 shows a sample frequency response plot for a natural frequency of 0.5 Hz and several damping ratios. 

With proper selection of the damping ratio, this transfer function allows the signal at low frequencies to be passed 
essentially unaltered. At the same time, the signal at frequencies beyond 0) n is attenuated greatly with a magnitude 
attenuation approaching 40 dB/frequency decade. Because equation (D-l) introduces a negative phase angle, it will 
have the effect of time-delaying the output signal. This lag is given by 

X = lim (D-4) 

to 0 CO CD n 



Figure D-2. Frequency response of low-pass infinite impulse response filter with time shifting for several damping 
ratios, w n = 0.5 Hz. 
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Because filtered output signals are to be time-correlated with unfiltered input signals, this time delay is accounted 
for by after-the-fact time-shifting of all filtered signals. This time-shifting is taken into account in figure D-2 which 
shows nearly zero phase lag for frequencies less than 0)„ . 

DIFFERENTIATOR TRANSFER FUNCTION 

The differentiator is derived from equation (D-l) by simply concatenating the filter transfer function with an 
open-loop zero (multiplying by s), which from Laplace transform mapping rules has the effect of differentiating the 

original signal. That is. 


= sl(s)-i( 0) 
_at 


The resulting transfer function is as follows: 


0(s) _ 

I(s) 



(D-5) 


(D-6) 


while using i(0) = 0 . Because the open-loop zero tends to infinite magnitude at high frequency, it must be 
concatenated with the low-pass filter to attenuate high-frequency measurement noise. Failure to do so results in a 
differentiated signal that is overwhelmed by the overamplified noise. 

Figure D-3 shows a sample frequency response plot for the differentiating filter at a natural frequency of 0.25 Hz 
and a damping ratio of 0.7071 with time-shifting. Notice that the magnitude appears to follow the proper slope 
through approximately 0.2 Hz. Beyond this frequency, the value rolls off, and the data no longer accurately 
represent the derivative. As a consequence, the derivative signal will always be reduced in frequency content from 
the original signal. As with the filtered data, a time delay occurs. This time delay is accounted for by after-the-fact 
tune-shifting of all differentiated signals, as shown by the nearly flat 90° phase at the lower frequencies. 


FILTER DISCRETIZATION 

The frequency-domain transfer functions are for continuous-time signals. For application to discrete-time- 
sampled data signals, the filters must be mapped to the discrete-time Fourier plane, z-plane, and inverse 
transformed to give difference equations. These difference equations can be implemented recursively to process the 
input signals. 

Map ping from the continuous-time frequency plane to the discrete-time Fourier plane is completed through the 
bilinear transform 6 . This mapping function is given by 


s = 



(D-7) 
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Figure D-3. Frequency response of low-pass infinite impulse response differentiator filter with time shifting, 
w„ = 0.25 Hz. 


where At is the sample interval of the discrete-time system, and z is the discrete-time Fourier transform variable. If 
z is substituted into equation (D-l) and the results are collected, the resulting z-transform function for the low-pass 
filter is as follows: 


0(z) _ b o(z+l) 2 

I(z) ~ 2 

v ' c "' a 2 z + 0/2 + 00 


(D-8) 
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where 


b 0 = (co„AO 

a 0 = 4 + (a) n Ar) 2 -4^co n Ar 
a x = - 8 + 2(to„A0 2 
a 2 = 4 + (co„A/) 2 + 4^co n A/ 


(D-9) 


Applying the rules of the inverse z-transform, the associated time-domain difference equation is as follows: 


°k = 


a 2 


(D-10) 


Similarly, if z is substituted into equation (D-6), and the results are collected, then the resulting z-transform 
function for the differentiating filter is as follows: 


O(z) _ c o(* - 1 ) 

a 2 z 2 + a x z + a 0 


(D-l 1) 


where 


c 0 = 2AtG> n 2 


(D-l 2) 


Applying the rules of the inverse z-transform, the associated time-domain difference equation is 


do\ 

dt) k 


r • • i ( d°\ ( do\ 


(D-13) 


F/ piatinns (D-10) and (D-13) are time-recursive equations with k being the discrete-time index. Because the 
equations are second-order, results from the two previous recursions are required at each frame. For initial startup, 
this condition will cause brief transients which damp out in a period of time on the order of several x. 

SECOND-ORDER ACCURATE, BACKWARDS-DIFFERENCE 
DIFFERENTIATOR 


If no filtering of the velocity or accelerations is desired, the HR differentiator described in the Differentiator 
Tr ansf er Function and Filter Discretization sections will not work. In this case, a second-order accurate, 
backwards-difference differentiator with the following form is used: 7 


do\ _ 3*‘jfc~4 ij^_i + i k _ 2 
dt J k 2 At 


(D-14) 
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VELOCITY AND ACCELERATION COMPONENTS 


The range, azimuth, and elevation angles may be filtered in the program. The geocentric position components of 
the vehicle, xc, yc, and zc, are differentiated and filtered once for velocity components and once again for 
acceleration components. The velocity and acceleration components are then rotated into the local north, east, and 
down system of the vehicle by 


- — 


- - 



v x 

K 

= r 

V y 

Yd 


1 


(D-15) 


Acc n 


Acc x 

Acc e 

= r 

ACCy 

Y cc d_ 


Acc z 


(D-16) 


where 


r = 


-sin(X,)cos(0) -sin(X,)sin(0) cos(A,) 
-sin(0) cos(0) 0 

-cos(X)cos(0) -cos (A,) sin (0) -sin(X) 


(D-17) 


and where X and 0 is the geodetic latitude and longitude of the vehicle, respectively. If the user desires, acceleration 
of gravity at the vehicle altitude, g, may be subtracted from downward acceleration. The value of g can be 
calculated from 16 


R, 


8 = 8 0 - 


(R 0 + zgeoidY 


(D-18) 


where g a is standard acceleration of gravity, R a is a radius of the Earth, and zgeoid is the geoid altitude of the 
vehicle. This subtraction makes the acceleration similar to that value measured by an onboard accelerometer. 
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APPENDIX E 


EARTH RELATIVE AND AIRDATA PARAMETER DERIVATION 

The Earth-relative parameters are a function of the Earth-relative velocity components V n , V e , and V d . Total 
Earth-relative velocity, V, is given by 


V = Jv „ 2 + V e 2 + V} (E-l) 

Flightpath heading, 'F, is the horizontal angular direction the vehicle is going relative to true north. This heading 
has nothing to do with where the vehicle is pointed and is computed as follows: 


¥ = 



(E-2) 


Flightpath angle, y, or the vertical angular direction that the vehicle is going relative to the local horizontal is 
defined positive upward as 


7 = 


tan 


-l 


-V A 




(E-3) 


Airdata parameters are determined from Earth-relative parameters that are adjusted for atmospheric properties. 
These atmospheric properties have been entered into a table input to the “radar” program and are discussed in 
appendix F. For interpolation and computation purposes, the magnitude and direction of the wind and lateral pres- 
sure altitude gradient are converted into Cartesian coordinates by 

W„ = IVcosCFj 
W t = W'sin('FJ 

= AH p i ar cos('¥ p i aI ) (E-4) 

= A H p ; a ,sin(4 , p ,„) 


The true airspeed, , of the vehicle is the square root of the squares of Earth relative speed plus the windspeed 
expressed in feet per second for each component. Vertical winds are assumed to be zero. 


V- = ■l(V„ + W„) 2 + (.V,+W'f + V/ 


(E-5) 
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True Mach number, , is the true airspeed divided by the speed of sound. That is. 


M m = 


k ]000R( T ~ + 273 i 5 o ) 


1200 m 

3937 ft 


(E-6) 


where T m is the ambient temperature at the altitude of the vehicle; k is the ratio of specific heats for air; R is the 
universal gas constant; and M a is the molecular weight of dry air. Pressure altitude from the atmospheric table is 
adjusted for a lateral pressure gradient by 


(E-7) 


where xr and yr are determined from equations (A-25) and (A-26). 

Ambient pressure, Ps M , is determined from pressure altitude by the U. S. Standard Atmosphere. 16 For regions of 
the atmosphere where the temperature has a constant lapse rate, L, the ambient pressure is 


Ps„ = Ps„ 


base 


1 + 


( r ^ + 273 - 15 °) 


\ P Phase ) 


-goK 

RL 


(E-8) 


where the subscript base denotes values at the lowest altitude of the region. For regions with a constant tempera- 
ture, ambient pressure is 


Ps„ = Ps _ 


( H u 'll 

*iT_ + m.\y ) \ n p- H p ba ")\ 


ban 


(E-9) 


The values of L, 7’ 00 ^^, and are given for each atmospheric region in table E-l. For altitudes greater than 

those in table E-l, ambient pressure is set to zero. 

Lastly, dynamic pressure, q , is defined as 


q = 


^MJ 


(E-10) 


43 



The values for wind magnitude, wind direction, and ambient pressure at the vehicle are determined by interpolation 
of the atmospheric table. 


Table E-l . Constants used in equations (E-8) and (E-9). 16 


Up, km 

L, K/km 

t„ ,°c 

°^base 

Ps ,lbf/ft2 

base 

0-11 

-6.5 

15.0 

2116.22 

11-20 

0.0 

-56.5 

472.68 

20-32 

1.0 

-56.5 

114.3450 

32-47 

2.8 

-44.5 

18.12885 

47-51 

0.0 

-2.5 

2.31630 

51-71 

-2.8 

-2.5 

1.398035 

71-84.852 

-2.0 

-58.5 

0.082631 
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APPENDIX F 


PROGRAM USE 

This appendix gives a short history of the “radar” program and operational details of its use. The primary focus of 
this discussion involves the program’s use at NASA Dryden. 

BACKGROUND 

An early version of this computer program dates back to the late 1950’s and was used with the X-15 program.* 
This code, version 1 .0, computed Cartesian position, altitude, total velocity, rate of climb, and radial acceleration. It 
used only the White Sands method for refraction correction (appendix C). Altitude was computed assuming a 
spherical Earth. Over the years, several new features and options were added to the program. In addition, the code 
was converted several times as mainframe computers were replaced with newer models. The program on the IBM 
computer was called version 2.0. Version 3.0 was on the CYBER computer. Version 3.0 was converted for use on 
the ELXSI computer as version 4.0. Version 5.0 was basically the current program implemented on the ELXSI, and 
version 6.0 was converted to run on a SUN 600 workstation. Versions 6.1 and 6.2 removed a few problems from 
version 6.0 and added a few new options. Previously, several versions of the same program existed simultaneously. 
In addition, each user tailored the code to fit individual needs. Inefficiencies and errors were found in the program, 
probably because of the lack of central control of the source code. To refurbish the program, each algorithm was 
investigated and derived from basic principles, and inefficiencies in the code were eliminated. Now the code, ver- 
sion 6.2, resides on the system level of csl, a SUN 600 workstation using System 4.1.2, and configuration changes 
will be tracked. 

EXECUTABLE AND SOURCE CODES 

For the computer csl at NASA Dryden, the files for the program reside in the directory, /dryden/cs 1/radar. The 
executable file radar handles input — output definitions and executes the FORTRAN code. The executable FOR- 
TRAN code, radar.source.6.2, is generated from the source code, radansource.6.2.f, by running radar.bind. 

PREPARING RAW DATA FOR THE SUN 600 COMPUTER 

The four possible sources for the raw radar data are as follows: 

• Flight Data Access System (FDAS) 2 at NASA Dryden. 

• Nine-track magnetic tape supplied by Western Aeronautical Test Range (WATR). 

• Engineering units range, azimuth, elevation in an unformatted file. 

• Geocentric Cartesian coordinate data in an unformatted file. 


*Marcalus, Vincent J., NASA-FRC Data Reduction Requirements for High Range Radar Data, undated working paper. 
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Figure F-l shows a flowchart for the radar program. The following subsections describe each source and provide 
specific directions for preparing the data for use on the SUN 600 computer called csl at NASA Dryden. 


■ Mandatory input and output files 
fii Optional 


Figure F-l. Flowchart of the “radar” program, version 6.2. 
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Flight Data Access System Data 

If the raw data resides on the FDAS, the user needs to copy the data into a working directory using the GetFdas 
utility. The steps to retrieve data using GetFdas are outlined next. 

1. Enter the utility with the command 


getfdas 


2. Select the flight project with the command 

project <flight-project> 

3. Enter the appropriate identification or password if one is required for your project. 

4. Select the appropriate flight with the command 

flight <flight-number> 

5. Complete the parameters command with 

parameters range_cts_## azimuth_cts_## elevation_cts_## 


where ## is the radar site number. Some vehicles are tracked by multiple radars simultaneously, and some projects 
have the main vehicle tracked by one radar and the chase vehicle tracked by a different radar. The “radar” program 
can analyze the data set from only one radar at a time. Valid radar site numbers for Edwards AFB, California, 
are 34 (primary NASA Dry den radar), 05 (backup NASA Dryden radar), and 38 and 41 (Air Force Flight Test Cen- 
ter, radars). 

6. Enter the write command with 


write <prefix>.raw.radar unc3 

where <prefix> is a descriptive name for the aircraft, flight, and maneuver, for example, fl04.1228.ad. This com- 
mand produces the file <prefix>.raw.radar in your working directory in unc3 format. 2 

7. Enter the time command with 

times hh.mm.ss.msec - hh.mm.ss.msec sync=radar_## 

where the first hh.mm.ss.msec is the hours, minutes, seconds, and milliseconds of the beginning of the data of 
interest, and the second time is the stop time. Some of the data on the FDAS are in Greenwich mean time (G.m.t.), 

♦Maine, Richard E., User’s Manual for GetFdas; Version 0.72, Apr. 30, 1993, working paper. 
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but the majority of these data are in local times. A maximum of 75,000 time points may be used for each run. This 
maximum corresponds to 60 min of data at 20 samples/sec with extra padding for filtering start-up transients, lags, 
and spike removal. As stated in appendix D, start-up transients die out on the order of several lag constants. Spike 
removal and filtering any parameter causes samples to be lost at the end of the data. Such losses result from the 
time shifting of the lag (equation (D-4)); therefore, the requested start time of the data should be before die require 
start time. In turn, the stop time of the data should be after the required stop time when filtering is performed. For 
the default-filtering frequencies shown in table F-l, 3.150 sec of extra data are needed to account for the lags. Spike 
removal requires one-half of a window of extra data points, which is 2.500 sec for the default window size and 

input sample rate. 


8. Exit the GetFdas utility and return to csl’s operating system by entering the command 


quit 


Table F- 1 . Namelist parameters for <prefix>.radansetup file. 


Namelist 

Parameter 

Description 

Type 

Default 

date 

month 

Month of flight 

Int*4 

None 


day 

Date of flight 

Int*4 

None 


year 

Year of flight 

Int*4 

None 

input 

prefix 

File name prefix for all input-output files 

Char*66 

None 


is tart 

Start hr, min, sec, msec time 

Int*4(4) 

None 


istop 

Stop hr, min, sec, msec time 

Int*4(4) 

None 


izulu 

Number of offset hours, 7 for G.m.t to P.d.t., 8 for 
G.m.t. to P.s.t. 

Int*4 

None 

indat 

spikes 

.true, to remove spikes 

Logic*4 

.false. 


window 

Number of points in spike window 

IntM 

100 


sigma 

Number of standard deviations for spike removal 

Real*8 

3. 


hlv 

Hold-last-value for missing data 

Logic*4 

.false. 


xi 

Damping ratio of filters 

Real*8 

Jl/2 


wbl 

Position-filtering break frequency, Hz 

Real* 8 

0.5 


wb2 

Velocity-filtering break frequency. Hz 

Real*8 

0.25 


wb3 

Acceleration-filtering break frequency, Hz 

Real* 8 

0.125 


gravity 

Is g to be subtracted from acceleration down? 

Logic*4 

.true. 


spsin 

Sample rate of input 

Real*8 

20. 
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Table F-l. Concluded 


Namelist 

Parameter 

Description 

Type 

Default 


nprint 

Print-thinning factor 

Int*4 

40 

amb 

corref 

.true, for refraction corrections 

Logic *4 

.true. 


emin 

el below which gradient refraction is used 

Real *8 

7° 


tdry 

Radar site dry bulb temperature 

Real*8 

59 °F 


twet 

Radar site wet bulb temperature 

Real*8 

59 °F 


pamb 

Radar site atmospheric pressure 

Real* 8 

27.25 in. Hg. 


reft 

.true, if refraction table is used 

Logic *4 

.false. 


nref 

Number of entries in refraction table 

Int*4 

0 


zmin 

Altitude below which gradient reflection is not 
truncated 

Real*8 

10,000. ft 


Is 

Length of gradient refraction segments 

Int*4 

1,000 ft 


grellip 

.true, if the full ellipsoid Earth model is to be used 
in gradient refraction corrections 

Logic *4 

.false. 

radsite 

sitlat 

Radar site geodetic latitude 

Real*8 

34.96081° 


siting 

Radar site longitude 

Real*8 

-117.91150° 


sith 

Radar site ellipsoid altitude 

Real*8 

2563.200 ft 


sitgs 

Radar site geodetic separation 

Real*8 

-99.393 ft 


zbias 

Altitude bias, zgeoid = h - sitgs - zbias 

Real*8 

0 


a 

Semimajor axis of Earth 

Real*8 

20925604.47 ft 


b 

Semiminor axis of Earth 

Real*8 

20855444.88 ft 


mlas 

Maximum mislevel, arc sec 

Real*8 

0 


mlsir 

Azimuth of maximum up mislevel, deg 

Real *8 

0° 

opt 

atm 

.true, if atmospheric table used 

Logic*4 

.false. 


numbp 

Number of entries in atmospheric table 

Int*4 

0 


binout 

.true, if Real*8 binary output desired 

Logic*4 

.false. 


binraw 

.true, if Real*8 binary raw t, r, az, el data are to be 
read in 

Logic *4 

.false. 


tape raw 

.true, if data on WATR 9-track tape 

Logic*4 

.false. 


unc3raw 

.true, if data on FDAS, unc3 format 

Logic*4 

.true. 


xyz 

.true, if Real *8 binary raw t, xc, yc, zc data are to 
be read in 

Logic *4 

.false. 


thin 

Thinning factor for output 

Int*4 

1 
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Western Aeronautical Test Range Nine-Track Tape Data 

If the radar data are from Edwards AFB, California, but are not on the FDAS, then the raw radar data are proba- 
bly on 9-track magnetic tape in encoded binary format supplied by the WATR. Most of these topes were written at 
a density of 800 bpi, and the rest were written at a density of 6250 bpi. The tape drive on the computer named cs2 
will automatically select the proper tope density. These data need to be copied from magnetic tope to the working 
directory. Steps for copying these data are as follows. 


1. Call the computer operator. 

2. Ask to have the tope mounted on the cs2 tope drive. 

3. Give the operator the tope number, for example, 3679f. 

4. Be sure to specify “no ring” to make the tope write protected. 

5. Read the tope using the command 

rsh -I tapeuser -n cs2 dd if=/dev/rst3 ibs=900 > <prefix>.raw.radar 


from the working directory. 

6. Verify that this procedure produced the file <prefix>.raw.radar in the working directory. 

7. Call the operator back to have the tape removed from the tope drive. 


Raw Range, Azimuth, Elevation Data 

If the data come from a radar that is not located at Edwards AFB, California, provide a binary file with time (sec- 
onds after midnight), range (feet), azimuth (degrees clockwise from true north), and elevation (degrees above local 
horizon) on each record and without a header. These four parameters must be in 8-byte unformatted words, as is 
generated by an unformatted write command in FORTRAN with double precision variables. Name this file <pre- 
fix>.raw.radar located in the working directory. 

Raw Cartesian Coordinate Data 

If the raw data are Cartesian coordinate position data referenced to the center of the Earth, provide an 8-byte 
unformatted word file with time (seconds after midnight), xc, yc , and zc (feet) on each record with no header. This 
coordinate system is described in appendix A. Name this file <prefix>.raw.radar located in the working directory. 

As shown in figure F-l, one of the four input sources named <prefix>.raw.radar is required. The file 
<filename>.radar.setup controls which of the four sources is used for a particular run. 

CHOOSING RADAR PROGRAM OPTIONS 

Generate the <prefix>.radar.setup file which contains all the option settings to execute the radar program using 
namelist and tobies. An example file is given at the bottom of this subsection. Ensure that the first column is blank. 
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The first record of this file is an 80-character line of text that will be the header for printed output. The namelists 
start on the second record. Table F-l lists these namelists, parameters, descriptions, and default values. Note that 
namelists “date” and “inpt” have no defaults and must always be included in the file. The parameter prefix should 
be set to <prefix>. This setting allows the program to find the appropriate input files. The refraction table, if one is 
desired (that is, if reft = .true.), follows the namelists. If airdata parameters are desired (that is, if atm = .true.) the 
atmospheric table is last in the file (fig. F-l). The namelist parameters and the table formats are discussed next. 

“date”, ‘inpt”, and ‘indat ” — The parameters of these namelists are adequately listed in table F-l, but a few 
caveats exist for the times and filtering break frequencies. There is a maximum of 75,000 time points/run, which 
corresponds to 60 min of data at 20 samples/sec with extra padding for filtering. If an interval is requested that 
exceeds the 75,000 time points, an error message will be printed and the program will stop. The start time should 
take into account start-up transients of the filters. The stop time is the actual desired stop time. The program will 
automatically read in enough data after the stop time to account for the filtering lags. 

“izulu” — If the raw data are from the WATR 9-track tape, time is in G.m.t., and izulu needs to be set to the hour 
offset between G.m.t. and the desired time zone. Data on the FDAS are generally in local time, so izulu would be 
zero. 

hlv — If gaps exist in the input data, the program will fill in the missing data. Data are generated by using hold- 
last-rate on the previous data as the default. The missing data can be generated by hold-last-value by setting hlv 
true, but remember that the velocities and accelerations will be adversely affected. 

wbl — If the position filtering break frequency, wbl, is set to zero, no filtering is performed. 

wb2 or wb3— If the break frequencies for either velocity or acceleration, wb2 or wb3, are set to zero, a 
second-order accurate, backwards-difference differentiator is used with no filtering for that derivative. 

corref — The parameter corref determines if atmospheric refraction corrections (both gradient refraction and 
White Sands methods) are performed. 

emin — Parameter emin is the measured elevation angle below which the gradient refraction method is used. 

If emin = -90°, use White Sands method. 

If emin = 90°, use the gradient refraction method. 

Giving an intermediate value to emin causes the White Sands method to be used for high-elevation angles where it 
is quite accurate, and the gradient refraction method to be used for low-elevation angles. Using the New Edwards 
curve, figure 4 can be used for Edwards AFB, California, to decide what value of emin is appropriate for a particu- 
lar application. Note that errors will be larger than those in figure 4 when the range exceeds 600,000 ft. If refraction 
corrections are to be performed, either the radar site atmospheric properties or a refraction table must be provided. 
If both are in the <prefix>.radar.setup file, the refraction table is used. 

tdry, twet, and pamb— The radar site atmospheric parameters are tdry, twet, and pamb. If the values of these 
three parameters are physically impossible, a warning will be printed, and the program will stop. 
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reft and nref— Parameters reft and nref are set if a refraction table is included for input. This table will have a 
column of geometric geoid altitude (above mean sea level) and a second column of refractivity (sometimes called 
“N” units). Table F-2 gives an example of a refractivity table. This table must be in increasing order of altitude and 
follow the namelists (fig. F-l). Because the program extrapolates refractivity when the desired value is outside the 
bounds of the table, the first two and last two entries in the table must exhibit the same trend as the rest of the data. 


Table F-2. Example of refractivity table (Pt. Arguello, California, 
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This information is usually given in 1000-ft increments. The White Sands method only uses the refractivity at the 
altitude of the radar site. 


zmin. Is, and grellip— The parameters zmin, Is, and grellip are used to control the gradient refraction method. 

• The zmin is the geoid altitude below which the truncation algorithm is not used. This truncation algorithm 
assumes an exponential decay of refractivity with altitude. Near the ground and especially in inversion layers, 
the refractivity may be irregular (fig. 2). Appendix C discusses this subject in additional detail. 

• The Is is the approximate length of the segment used in the gradient refraction method. The accuracy of the 
gradient refraction method depends on Is, and its optimum value depends on the roundoff and truncation errors 
of the computer. For the SUN 600 computer at NASA Dryden, a value of Is = 1000 ft is generally used although 
this value may not be the optimum. This value was selected because data from balloon profiles of refractivity 
come at 1000-ft intervals. 

• The grellip controls the Earth model used during the gradient refraction method. If grellip = .false., the Earth 
model is spherical. If grellip = .true., the Earth model is the ellipsoid. To date, no significant differences in 
refraction corrections using the two settings of grellip have been observed, but possibly it could be significant 
at extreme ranges. Using the ellipsoid model during the gradient refraction method increases computation time 
tremendously. 

sitlat, siting, sith, and sitgs — These parameters represent the geodetic latitude, longitude, ellipsoid altitude, and 
geodetic separation, respectively, of the radar site for the Earth model described by the parameters a and b. The 
defaults are for radar 34 at Edwards AFB, California, in the WGS 84 system. Table F-3 gives radar site positions for 
several installations in the WGS 84 system. Table 1 gives values for a and b in several systems. Entering values for 
sitlat, siting, sith, sitgs, a, or b has no effect for the FDAS or WATR input formats. When the raw data are from the 
FDAS or the WATR 9-track tape, the radar site coordinates in the WGS 84 system are automatically entered to the 
program for any of the four radars at Edwards AFB, California. 


Tabl^F^^Califomia radar site coordinates in WGS84 system. 


Radar site 

Latitude, 

deg 

Longitude, deg 

h s ,ft 

gs, ft 

zgeoid s , ft 

Edwards 34 

34.96081 

-117.91150 

2563.200 

-99.393 

2662.593 

Edwards 05 

34.95774 

-117.91187 

2538.699 

-99.403 

2638.102 

Edwards 38 

34.96961 

-117.92941 

2605.359 

-99.304 

2704.663 

Edwards 41 

34.97045 

-117.93056 

2623.771 

-99.298 

2723.069 

Pt. Mugu 003004 

34.12289 

-119.15475 

-77.211 

-121.027 

43.816 

Pt. Pillar 213002 

37.49784 

-122.49970 

61.030 

-108.917 

169.947 

Pt. Pillar 213003 

37.49687 

-122.49667 

7.192 

-108.881 

116.073 

San Nicholas Is. 013003 

33.24769 

-119.52074 

809.700 

-122.040 

931.740 

Vandenberg 023001 

34.58276 

-120.56157 

2059.281 

-113.510 

2172.791 

Vandenberg 023002 

34.58305 

-120.56111 

2059.310 

-113.500 

2172.811 

Vandenberg 023003 

34.66586 

-120.58144 

288.871 

-113.412 

402.283 

Vandenberg 033001 

34.75825 

-120.62712 

85.298 

-113.103 

264.019 

Vandenberg 033701 

34.77492 

-120.53607 

404.530 

-112.674 

517.204 
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zbias— The user can bias the altitude using the zbias parameter. This function may be desirable if the geoid sep- 
aration at the radar site differs significantly from the geoid separation at the location of the vehicle. To output alti- 
tude as ellipsoid altitude, make zbias = -sitgs. 

an d mldir — These parameters allow corrections for radar pedestal mislevel. It is assumed that the pedestal 
is tilted by the number of arc seconds in mlas and that the azimuth of the high side of the pedestal is given by 

mldir. 


atm This parameter is set to .true, if an atmospheric table is input so that airdata parameters can be calculated. 

This table consists of seven columns, an example of which is given in table F-4, and is placed at the end of the file 
(fig. F-l). The first column is geometric altitude (above mean sea level) in feet, zgeoid, and the second is the corre- 
sponding pressure altitude in feet, H p ta ^i e > for a weather analysis above the radar site. Ambient temperature, T x , 
is shown the third column. The fourth and fifth columns include windspeed in knots, IV, and wind direction, 'F w • 
The wind direction is in degrees from true north and describes the direction from which the wind blows, so a north 
wind blows from north to south, and y ¥ w = 0°. The sixth and seventh columns list lateral pressure altitude gradi- 
ent and direction, A H lat and lal . The lateral pressure altitude gradient is the number of feet that pressure 
altitude changes per nautical mile, and direction is the direction of decrease of pressure. As a result, if a low pres- 
sure area is due north of the radar site, then V p lat = 0 C . This table must be arranged in increasing order of geo- 
metric altitude. If values outside the range of the table are sought, the nearest values of the table are used. 

numbp— This parameter represents the number of entries (rows) in the atmospheric table. 


Table F-4. Sample atmospheric table for radar 34, at Edwards, California, on June 13, 1988. 

— — — — — — e ^»aeg= 







AH p lat ’ de S 

Vp lat ’ deg 

zgeoid , ft 

H p table * ft 

U 

# . 

h, 8 

W.kn 

4V deg 

ft/n. mi. 

2,302 

2,275 

30.0 

03 

085 

0.10 

250 

3,330 

3,243 

25.5 

04 

074 

0.10 

170 

4,133 

4,003 

22.7 

10 

074 

0.10 

165 

4,956 

4,781 

21.5 

13 

080 

0.20 

160 

6,145 

5,902 

20.0 

12 

072 

0.20 

160 

10,354 

9,882 

8.8 

09 

145 

0.25 

235 

13,597 

12,977 

-0.6 

06 

275 

0.25 

360 

15,806 

15,079 

-0.2 

11 

310 

0.30 

040 

19,201 

18,289 

-8.1 

14 

313 

0.40 

057 

24,759 

23,574 

-21.4 

20 

306 

0.50 

032 

28,800 

27,469 

-31.5 

43 

300 

0.55 

032 

31,495 

30,065 

-38.3 

34 

300 

0.55 

032 

33,045 

31,570 

-42.0 

33 

301 

0.65 

030 

35,533 

33,999 

-48.8 

38 

300 

0.67 

027 

38,258 

36,679 

-55.0 

39 

300 

0.70 

020 

39,232 

37,647 

-56.5 

44 

294 

0.74 

020 

40,244 

38,662 

-58.8 

50 

300 

0.76 

020 

41,855 

40,284 

-59.4 

45 

301 

0.74 

020 
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The program uses the distance north and east of the radar site, xr and yr, from equations (A-25) and (A-26). 
Because these quantities originate at the radar site, the atmospheric table should be referenced to the specific radar 
site geographically. If the atmospheric table is not referenced to the radar site, the radar data must first be processed 
without atmospheric input to give geocentric position. Then using the geocentric position as input, run the program 
again with the radar site respecified as the origin of the atmospheric data. 

binraw, taperaw, unc3raw, xyz, and binout — These parameters are explained in table F-l and control which 
input source is used and whether binary output is generated. Note that although binary output takes up a great deal 
more memory than the normal compressed output (which is in cmp3 format) 2 , it retains its full 8-byte word. The 
cmp3 format output retains only the three most significant bytes. This level of precision is acceptable for most 
parameters. On the other hand for such parameters as Cartesian position from the center of the Earth, latitude, and 
longitude, a significant amount of resolution is lost when the compressed format is used. 

thin — This last parameter in the namelist is the requested thinning factor for the output. This parameter has no 
effect on the filters or differentiation and is only used for the output routines. 

Figure F-2 shows an example radar.setup file. This file will process data for flight 1228 of the F104 aircraft 
(Lockheed Corporation, Burbank, California) which occurred on June 13, 1988. Four and one-half minutes of data 
are analyzed from 1 1:33:32 to 1 1:38:02 local time, and a 7-hr offset exists between P.d.t. and G.m.t. The filtering 
frequencies have been changed from the defaults, and the radar site weather conditions are entered for atmospheric 
refraction corrections. Radar site information will be entered automatically when the raw data are read, but manual 
entry of the namelist name is still required in this setup file. An atmospheric table has been provided with 1 8 break- 
points, so airdata quantities will be computed. The input is from the WATR 9-track tape, and a binary output file 
will be generated. Entering the command 


man radar 

on csl gives an abbreviated set of these instructions and this example. 

RUNNING THE RADAR PROGRAM 

Now that the input files have been created, begin the radar program by typing 

radar <prefix>.radansetup 

in the working directory. If the program will take a while to run (for example, if the gradient refraction is used), 
then consider using this program in another window or in the background. After successful use of the program, the 
compressed output in cmp3 format 2 will be in <prefix>.radar,out. The binary output, if requested, will be in 
<prefix>.radar.out.bin (fig. F-l). If an atmospheric reference table is used, all parameters given in table F-5 are 
output. If an atmospheric reference table is not provided, only channels 0 through 29 are output. Most of the output 
parameter names have an “r” prefix to help distinguish these radar-derived results from results derived from 
aircraft-calibrated parameters, global-positioning systems, or inertial navigation systems. The data can be exam- 
ined using a plotting program. A text output will be created called <prefix>.radar.print. An abbreviated version of 
the text output is sent to standard output, which is the monitor when the program is run interactively. Multiple data 
sets may be analyzed simultaneously from the same directory if each has a unique prefix. 
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radar data for F104 fltl228 a/c, radar #34 

$date month=06 , day=13, year=1988 $ . __ , 

$inpt prefix=' f 104. 1228' , istart=ll , 33 , 32 , 000, istop=ll , 38 , 02 , 000 , izulu-7 $| 
$indat wbl=0.2, wb2=0.1, wb3=0.05 $ 

$amb tdry=86 . 0 , twet=59.0, pamb=27.17 $ 

$radsite $ 

$opt atm= . true . , numbp= 18 , binout= . true . , taperaw= . true . $ 


2302. , 

, 2275. 

, 30.0, 

03. , 

, 085. 

, 0.10, 250. 

3330. , 

, 3243. 

, 25.5, 

04. 

, 074. 

, 0.10, 170, 

4133. 

, 4003. 

, 22.7, 

10. 

, 074. 

,0.10, 165. 

4956. 

„ 4781. 

, 21.5, 

13. 

, 080. 

,0.20, 160 

6145. 

, 5902. 

, 20.0, 

12. 

, 072. 

,0.20, 160 

10354. 

, 9882. 

, 8.8, 

09. 

, 145. 

,0.25, 235 

13597. 

, 12977. 

, -0.6, 

06. 

, 275. 

, 0.25, 360 

15806. 

, 15079. 

, -0.2, 

11. 

, 310. 

,0.30, 040 

19201. 

, 18289. 

, -8.1, 

14. 

, 313. 

,0.40, 057 

24759. 

, 23574. 

, -21.4, 

20. 

, 306. 

,0.50, 032 

28800. 

, 27469. 

, -31.5, 

43. 

, 300. 

,0.55, 032 

31495. 

, 30065. 

,-38.3, 

34. 

, 300. 

,0.55, 032 

33045. 

, 31570. 

,-42.0, 

33 . 

, 301. 

,0.65, 030 

35533. 

, 33999. 

, -48.8, 

38. 

, 300. 

,0.67, 027 

38258. 

, 36679. 

, -55.0, 

39. 

, 300. 

,0.70, 020 

39232. 

, 37647. 

, -56.5, 

44. 

, 294. 

,0.74, 020 

40244. 

, 38662. 

, -58.8, 

50. 

, 300. 

,0.76, 020 

41855. 

, 40284. 

, -59.4, 

45. 

, 301. 

,0.74, 020 


Figure F-2. Sample <preflx>.radansetup file for flight 1228 of the F104 aircraft, June 13, 1988. (The first col- 
umn of this file is blank.) 


MEMORY 

On the SUN 600 computer, the FORTRAN source code takes up 89 kilobytes of memory, and the compiled code 
takes up 385 kilobytes of memory. For 30 min of data at 20 samples/sec, the approximate maximum file sizes are as 
follows: 

• 84 kilobytes for the standard output saved to a file 

• 656 kilobytes for <prefix>.radar.print 

• 3.2 megabytes for <prefix>.radar.out 

• 1 1 megabytes for <preflx>.radar.out.bin 

• 23 megabytes for the assorted scratch files. 

Farh scratch file is automatically deleted when the program finishes using that file. File sizes are smaller when 
shorter time segments are used and if the airdata parameters are not desired. The <prefix>.raw.radar file may 
require several megabytes. Table F-5 shows the output parameters. 
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Table F-5. Output parameters. 


No. 

Parameter 

Description 

Equation 

0. 

time 

Local time after midnight, sec 

— 

1. 

rraw 

Range, counts 

— 

2. 

araw 

Azimuth, counts 

— 

3. 

eraw 

Elevation, counts 

— 

4. 

reng 

Range, no corrections, ft 

— 

5. 

aeng 

Azimuth, no corrections, deg 

— 

6. 

eeng 

Elevation, no corrections, deg 

— 

7. 

rfilt 

Range, spikes removed and filtered, ft 

— 

8. 

afilt 

Azimuth, spikes removed and filtered, deg 

— 

9. 

efilt 

Elevation, spikes removed and filtered, deg 

— 

10. 

rcor 

Range, rfilt corrected for refraction, ft 

— 

11. 

ecor 

Elevation, efilt corrected for refraction, deg 

— 

12. 

rxc 

Distance from center of Earth along equitorial plane toward 0° longitude, 

ft 

(A-20) 

13. 

rye 

Distance from center of Earth along equitorial plane toward 90° east 

longitude, ft 

(A-20) 

14. 

rzc 

Distance from center of Earth toward north pole, ft 

(A-20) 

15. 

rxr 

Distance north of radar site on ellipsoid Earth, ft 

(A-24) 

16. 

ryr 

Distance east of radar site on ellipsoid Earth, ft 

(A-25) 

17. 

rz 

Geometric altitude above mean sea level (geoid), zgeoid = h- sitgs - 

zbias, ft 

(B-40) with bias 

18. 

rglat 

Geodetic latitude, deg 

(A-22) 

19. 

rgclat 

Geocentric latitude, deg 

(A-21) 

20. 

rglong 

Longitude, deg 

(A-23) 

21. 

rvn 

Velocity north filtered, ft/sec 

(D-15) 

22. 

rve 

Velocity east filtered, ft/sec 

(D-15) 

23. 

rvd 

Velocity down filtered, ft/sec 

(D-15) 

24. 

ran 

Acceleration north filtered, ft/sec 2 

(D-16) 

25. 

rae 

Acceleration east filtered, ft/sec 2 

(D-16) 

26. 

rad 

Acceleration down filtered, local gravity may be subtracted, ft/sec 2 

(D-16 and D-18) 

27. 

rvtot 

Earth relative total velocity, ft/sec 

(E-l) 

28. 

rfjph 

Flightpath heading, deg 

(E-2) 

29. 

rfpa 

Flightpath angle, deg 

(E-3) 

Atmospheric-derived parameters 

30. 

rwindmag 

Wind magnitude, ft/sec 

— 

31. 

rwinddir 

Wind direction from true north, deg 

— 

32. 

rpsinf 

Ambient pressure, lbf/ft 

(E-8 and E-9) 

33. 

rtinf 

Ambient temperature, °R 

— 

34. 

rvinf 

True airspeed, ft/sec 

(E-5) 

35. 

rminf 

Mach number 

(E-6) 

36. 

rqbar 

Dynamic pressure, lbf/ft 

(E-10) 

37. 

rhp 

Pressure altitude, corrected for lateral pressure gradient, ft 

(E-7) 
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APPENDIX G 


NOMENCLATURE 


A 

AFB 

AFFTC 

AN/FPS-16 

Acc 


a 

a, b, c, d,- f, g 


az 

a 0 , flj, fl2> K c 0 

a2, b2, c2 



C c» d c 


bpi 


b q , c q , d q , e q 


C 

C 

Char 


D 

d 

delt 


e 

el 

es 

ev 

FDAS 

G.m.t. 

GPS 

8 


scale height constant for equation (C-7), m 
Air Force Base 

Air Force Flight Test Center, Edwards, California 

Army-Navy/Fixed Position System radar 

vehicle acceleration, ft/sec 2 

ellipsoid semimajor axis, ft 

refractivity constants for equations (C-3) and (C-4) 

azimuth angle, deg 

filter constants 

coefficients of second-order form of quartic equation 
scale height constant for equation (C-7), m 
ellipsoid semiminor axis, ft 
coefficients of resolvent cubic equation 
bits per inch 

coefficients of quartic equation 
dummy variable 

scale height constant for equation (C-7), m 

character FORTRAN variable 

speed of light in a vacuum, ft/sec 

downrange distance from radar site, ft or yd 

distance from ellipsoid, ft 

time interval for gradient refraction segment, sec 

eccentricity 

elevation angle, deg 

saturation vapor pressure, in. Hg. 

vapor pressure, in. Hg. 

Flight Data Access System 

Greenwich mean time 

Global Positioning System 

acceleration of gravity at the vehicle, ft/sec 2 
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So 

8S 

H 

Hp 

Hp table 

h 

I 

HR 

INS 

Int 

i 

ia 

J 

j 

K, , K 2 , Kj , K 2 

l * x e l r 

k 

L 

L 

LSB 

Logic 

Ls 

Is 

M 

K 

Ms 

M n 

my n 

m.s.l. 

N 

NASA 

n 

ns 


standard acceleration of gravity, 32.1740 ft/sec 2 or 9.80665 m/sec 2 
geoid separation, ft 
scale height, m 
pressure altitude, ft 

pressure altitude from atmospheric table, ft 

geometric ellipsoid altitude, ft 

transformed input function 

infinite impulse response 

inertial navigation system 

integer FORTRAN variable 

time domain input function 

internal Earth angle, deg 

equation (B-2) 

imaginary number, J-\ 

White Sands method coefficients 
ratio of specific heats for air, 1 .4 
Laplace transform 
temperature lapse rate, K/km 
least significant bit 
logical FORTRAN variable 
lower segment length, ft 

gradient refraction segment length when T| = 1, ft 
magnitude 

molecular weight of dry air, 28.9644 gm/mole 

middle segment length, ft 

true firee-stream Mach number 

average remaining segment length, ft 

coefficients of linear solution to quartic equation 

mean sea level 

refractivity 

National Aeronautics and Space Administration 
segment number 
number of segments 



o 


o 

P.d.t. 

P.s.t. 

Ps„ 

ps 

P, Q 

q 

R 

Real 

Re 

*o 

r 

r 

rh 

S 


s 

s 

T, 

L 

Twet s 

t 

*c 

Us 

u 

V 

Voo 

v 


V 

W 

WATR 

WGS 


w 


transformed output function 

time domain output function 

Pacific daylight time 

Pacific standard time 

ambient pressure, lbf/ft 2 

radar site atmospheric pressure, in. Hg. 

coefficients of reduced cubic equation 

dynamic pressure, lbf/ft 2 

universal gas constant, 8.31432 kg-m 2 /(K-mole-sec 2 ) 
real FORTRAN variable 

local radius of curvature of the Earth at the radar site, ft 

radius of the Earth for g calculation, 20,855,531.5 ft 

ratio of reduced cubic and trigonometric reduced cubic equation solutions 

range, ft 

relative humidity, percent 
dummy variable 
radar site vertical vector 
Laplace transform variable 
dry bulb temperature, °R or °F 
ambient temperature, °C or °R 
wet bulb temperature, °R or °F 
time, sec 

variable of trigonometric reduced cubic equation 

upper segment length, ft 

ellipse major coordinate, ft 

Earth relative velocity, ft/s 

true airspeed, ft/sec 

local vertical vector 

ellipse minor coordinate, ft 

windspeed, ft/sec or kn 

Western Aeronautical Test Range 

World Geodetic Survey 

width of radar segment (separation between beams), ft 
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X, y,z 
xc, yc, zc 

*t> y* *i 

x f 

xr 

yr 

z 

z 

ze 

zgeoid 

Subscripts 

base 

d 

e 

k 

m 

min 

n 

r 

s 

x, y, z 

Superscripts 


Symbols 

a 

p 

r 

Y 

A 


point on ellipsoid, ft 

Cartesian distance from the center of the Earth, ft 
radar site local north, east, and down coordinates, ft 
variable of quartic equation 

distance from the radar site north to the latitude of the piercing point, ft 
variable of resolvent cubic equation 

distance from the radar site east to the longitude of the piercing point, ft 

vertical distance from radar site, ft or yd 

discrete-time Fourier transform variable 

variable of reduced cubic equation 

geoid altitude estimate, ft 

geometric altitude above mean sea level (geoid), ft 

value at the lowest altitude of the region 

down component 

east component 

sample number 

measured 

minimum 

north component 

corrected for atmospheric refraction 
radar site 

geocentric coordinates 

piercing point coordinates from center of Earth, ft 

values scaled by 1/a 
approximated 

Lagrange multiplier 

reduced latitude, deg 

rotational matrix 

flightpath angle, deg 

discriminant of reduced cubic equation 



AH 
A el 

Ah n 

Ar 

At 

6 

^"ren 

T| 

e 

x 


p lat 


p lat 


% 

a 

T 

<t> 

4* 

4* 

4* 

w 

(0 

© 

z 

Namelists 

date 


day 

month 

year 

inpt 

is tart 


istop 

izulu 

prefix 


lateral pressure altitude gradient magnitude, ft/n. mi 
elevation angle error, army mils 

estimate of altitude of the target above the current segment, ft 

range error, yd 

sample interval, sec 

turning angle, deg 

total rentmining turning, deg 

average remaining turning angle, deg 

index of refraction 

longitude, deg 

geodetic latitude, deg 

geocentric latitude, deg 

damping ratio 

standard deviation 

lag, sec 

phase angle, deg 

flightpath heading from true north, deg 

direction of lateral pressure gradient (towards pressure decrease), deg 

wind direction from true north (from which wind blows), deg 

frequency, Hz 

natural frequency, Hz 

ellipsoid piercing point 

trigonometric solution angle, deg 


day of flight 
month of flight 
year of flight 

data start time 

data stop time 

number of offset hours 

prefix of the names of the radar files 
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indat 


gravity 

hlv 

nprint 

sigma 

spikes 

spsin 

wbl 

wb2 

wb3 

window 

xi 

amb 

corref 

emin 

grellip 

Is 

nref 

pamb 

reft 

tdry 

twet 

zmin 

radsite 

a 

b 

mlas 

ml di r 

sitgs 

sith 

sitlat 

siting 

zbias 


logical for local gravity calculation 
logical for hold-last-value for missing data 
printing thinning factor 

number of standard deviations for spike removal 
logical for spike removal 
sample rate of input 
position-filtering break frequency, Hz 
velocity-filtering break frequency, Hz 
acceleration-filtering break frequency, Hz 
number of points in spike removal window 
filter-damping ratio 

logical for refraction corrections 

namelist parameter regulating refraction corrections, deg 

logical for Earth model used in refraction correction 

segment length for gradient refraction method, ft 

number of rows in refraction table 

ambient pressure at radar site, in. Hg. 

logical for refraction table 

dry bulb temperature at radar site, °F 

wet bulb temperature at radar site, °F 

namelist parameter regulating truncation algorithm, ft 

ellipsoid semimajor axis, ft 
ellipsoid semiminor axis, ft 
maximum mislevel, arc sec 
azimuth of maximum up mislevel, deg 
geoid separation at radar site, ft 
ellipsoid altitude of radar site, ft 
geodetic latitude of radar site, deg 
longitude of radar site, deg 
altitude bias, ft 
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opt 


atm 

logical for atmospheric table 

binout 

logical for 8-byte word binary output data 

binraw 

logical for 8-byte word binary raw input data 

numbp 

number of rows in atmospheric table 

taperaw 

logical for WATR 9-track tape raw input data 

thin 

output-thinning factor 

unc3raw 

logical for Flight Data Access System input data, unc3 format 

xyz 

logical for Cartesian input data 
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